F/8  20/* 


AD— A 108  123 
UNCLASSIFIED 


FOREION  TECHNOLOOY  DIV  iR I 8HT -PATTERSON  AFB  OH 
ACTA  NECHANlCA  SINICA  (SELECTED  ARTICLES) . (U) 

OCT  81 

FTP-ID  (RS)T-0707-81 _ NL_ 


0HO  FILE  COEV  ADA108123 


FTD-ID (RS)T-0707-8l 


EDITED  TRANSLATION 


FTD-ID(RS)T-0707~8l  29  October  1981 

MICROFICHE  NR:  ftd-81-C-000989 

ACTA  MECHANICA  SINICA  (Selected  Articles) 

English  pages:  56 

Source:  Acta  Mechanlca  Sinica,  Nr.  4,  1979, 
pp.  303-322,  344-352 

Country  of  origin:  China 
Translated  by:  SCITRAN 

F33657-78-D-0619 
Requester:  FTD/TQTA 

Approved  for  public  release;  distribution _ 

unlimited.  ft-' 


THIS  TRANSLATION  IS  A  RENDITION  OR  THE  ORIGI- 
NAL  FOREIGN  TEXT  WITHOUT  ANY  ANALYTICAL  OR 
editorial  COMMENT.  STATEMENTS  or  theories 
ADVOCATEDOR  IMPLIED  ARE  THOSE  OR  THE  SOURCE 
AND  DO  NOT  NECESSARILY  RERLECT  THE  POSITION 
OR  OPINION  OR  THE  FOREIGN  TECHNOLOGY  Dl- 
VISION. 


PREPARED  BY: 

TRANSLATION  DIVISION 
FOREIGN  TECHNOLOGY  DIVISION 
WP-AFB,  OHIO. 


FTD 


Date  29  Oct  19  81 


TABLE  OF  CONTENTS 


Generalized  Variational  Principles,  Variational 
Principles  and  Dual  Extremum  Principles  of  Gas- 
dynamics  for  the  Direct  Problem  of  Turbumachine 
Airfoil  Cascades  on  an  Arbitrary  Stream  Sheet  of 
Revolution,  by  Liu  Gao-Lian . 

Calculation  of  Unstarted-Choked  Regime  of  Supersonic 
and  Transonic  Compressor  Cascades,  by  Cheng  Jia-gang 

The  Numerical  Computation  for  Diffraction  and 
Reflection  of  a  Traveling  Shock  Wave  From  an  Open 
End  of  a  Shock-Tube,  by  Du  Xun,  Wei  Zhong-lei, 

Li  Wen-xuan,  Wang  Jia-jun . 


i 


GENERALIZED  VARIATIONAL  PRINCIPLES,  VARIATIONAL  PRINCIPLES  AND 
DUAL  EXTREMUM  PRINCIPLES  OP  GASDYNAMICS  FOR  THE  DIRECT  PROBLEM 
OP  TURBOMACHINE  AIRFOIL  CASCADES  ON  AN  ARBITRARY  STREAM  SHEET 
OP  REVOLUTION* 


Liu  Gao-Lian 

ABSTRACT 


The  present  paper  is  devoted  to  the  search  for  variational  principles  for  the  di¬ 
rect  aerodynamic  problem  of  channels  and  aerofoil  cascades  (with  prescribed  distribu¬ 
ted  mass  injection  or/and  suction  along  the  profile  contours  for  blade  cooling  or 
boundary  layer  control)  on  an  arbitrary  stream  sheet  of  revolution  and  aimed 
mainly  at  providing  a  rigious  and  sound  theoretical  basis  for  introducing  and  widely 
applying  the  finite  element  method  to  computational  aerodynamics  of  turbomachines. 
Three  generalized  variational  principles  (among  which  two  are  for  homoentropic 
flow,  and  the  other  for  nonhomoentropie  flow)  are  developed  first,  and  then  from 
them  three  families  of  variational  principles  are  derived  respectively  by  means  of 
successive  transformations.  Furthermore,  for  homoentropic  subsonic  flows  a  pair  of 
dual  extremum  principles  are  also  shown  to  exist 

Special  attention  is  payed  to  taking  full  advantage  of  2  effective  means:  “na¬ 
tural  boundary  conditions"  and  “artificial  boundaries",  in  order  to  greatly  simplify 
the  numerical  treatment  of  various  boundary  conditions  characterizing  the  complcxitg 
of  the  cascade  problem.  So  in  all  variational  principles  all  boundary  conditions 
have  been  converted  to  natural  ones. 

The  variational  principles  presented  here  can  be  used  to  advantage  as  basis  also 
for  variational-finite  difference  method  and  tor  various .  direct  variational  methods 
(e.g.  Rita’s  method,  Kantorovich’s  method  etc.). 


received  on  February  26,  1978. 

■^This  paper  has  been  presented  in  the  National  Experience  Ex 
change  Meeting  in  the  thermodynamic  calculation,  design  and 
experiment  of  the  aerodynamics  of  turbomachines. 
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This  paper  is  devoted  to  the  investigation  of  the  variational 
principles  for  the  direct  problem  of  channel  and  airfoil  cascade 
aerodynamics  on  an  arbitrary  stream  sheet  of  revolution  (with  gas 
injection  and/or  suction  along  the  profile  contour  of  the  blade. 
For  example,  boundary  layer  control,  gas  film  or  perspiration 
type  of  cooling,  etc.)  in  order  to  provide  a  more  rigorous  theo¬ 
retical  basis  for  the  introduction  and  expansion  of  the  finite 
element  method  in  the  field  of  turbomachine  aerodynamic  calcula¬ 
tion.  The  results  include:  1.  For  homoentropic  flow,  two  var¬ 
iational  principles  were  established,  and  then  from  them  two 
families  of  semi-generalized  variational  principles,  including  a 
pair  of  dual  extremum  principles  under  subsonic  condition  were 
derived;  2.  For  nonhomoentropic  flow,  a  generalized  variational 
principle  was  established  and  then  from  it  a  family  of  sub-gen¬ 
eralized  variational  principles  was  derived.  This  paper  placed 
special  emphasis  on  taking  full  advantage  of  the  natural  boundary 
conditions  and  artificial  boundaries  to  simplify  the  treatment  of 
various  boundary  conditions  of  the  complex  cascade  problem.  In 
addition  to  being  in  the  finite  element  method,  they  can  be  used 
to  create  advantageous  conditions  for  the  variational-finite  diff¬ 
erence  method  and  for  various  direct  variational  methods. 

The  use  of  the  variational  principle  to  describe  a  physical 
process  has  the  unique  advantage  of  simple  straight  forwardness. 

It  can  use  its  unique  natural  boundary  conditions  to  simplify 
(.Including  on  the  discontinued  surface)  the  treatment  of  boundary 
conditions.  It  has  always  been  the  basis  of  various  direct  var¬ 
iational  methods.  The  generalized  variation  principle  emerged  in 
the  fifties  [1].  It  created  conditions  to  simplify  the  selection 
of  coordinate  functions  as  well  as  to  increase  the  accuracy  in 
the  direct  variation  method.  It  has  been  primarily  used  in  solid 
state  mechnaics.  The  investigation  and  applications  of  the  var¬ 
iational  princlpl-  in  fluid  mechanics  have  been  rare  [3-6].  Its 
major  reasons  most  probably  are  the  complexity  of  the  boundary 
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conditions  (especially  in  the  case  of  cascades)  and  the  non¬ 
linearity  of  the  problem.  In  the  recent  decade,  the  development 
of  the  finite  element  method  brought  new  interest  to  the  class¬ 
ical  variational  principles.  The  variational  principle  is  one 
of  the  theoretical  foundations  of  the  finite  element  method. 

The  characteristics  of  the  finite  element  method  (e.g.,  arbitrary 
slicing  and  slice  implantation)  also  greatly  overcomes  the  diffi¬ 
culties  encountered  in  the  treatment  of  the  complicated  boundary 
conditions  using  the  direct  variational  method.  Thus,  it  elim¬ 
inated  the  disadvantage  that  such  a  method  was  only  suitable  for 
the  general  characteristics,  but  not  suitable  for  the  field  dis¬ 
tribution.  Therefore,  the  combination  of  the  two  methods  is 
really  complementary. 

This  paper  for  the  first  time  established  the  variational 
principles,  the  variational  principle  families  and  the  dual  extre¬ 
mum  principles  of  gas  dynamics  for  the  directproblem  of  cascades 
and  channels  of  turbomachines  on  an  arbitrary  stream  sheet  of  revo¬ 
lution.  Taking  into  account  the  needs  of  the  newly  developed  area 
(such  as  boundary  layer  control,  gas  film  or  perspiration  type  of 
cooling,  etc.),  we  have  included  the  conditions  of  direct  prob¬ 
lems  with  gas  injection  and  suction  on  the  surface  of  the  blade. 

To  further  simplify  the  treatment  of  the  boundary  conditions  ,  spe¬ 
cial  attention  was  paid  to  fully  take  advantage  of  the  unique 
actions  of  the  natural  boundary  conditions  and  the  artificial 
boundaries.  The  latter  enabled  us  to  reach  a  better  coordination 
on  the  intersection  of  each  element  in  the  finite  element  method. 

In  the  classical  direct  method,  it  was  possible  to  select  coordi¬ 
nation  functions  more  flexibly  by  area,  which  simultaneously 
assured  the  continuity  on  the  artificial  boundary  to  the  extent 
possible. 


I .  The  gas  dynamic  equations  of  cascades  on  an  arbitrary  stream 
sheet  of  revolution 


We  limited  our  investigation  to  the  condition  that  the  gen¬ 
eralized  enthalpy  H  (.also  known  as  the  stagnant  evolution  enthal¬ 
py)  of  the  incoming  flow  is  homogeneous  Ci.e.,  V*=  0).  Here  H  * 
o  o  * 

ith  w*V2)-(u  /2)  and  V  represents  the  gradient  along  the  surface 

of  the  flow.  Therefore,  for  a  system  rotating  at  a  constant  ang- 
ular  velocity  yr,  the  gas  dynamic  equations  [ 7 J  for  an  ideal  gas 
flowing  along  an  arbitrary  stream  sheet  of  revolution  under  steady 
and  adiabatic  conditions  can  be  written  in  the  following  dimen¬ 
sionless  forms  (Figure  1): 

1 .  Homoentroplc  flow  (Figure  2(b) 

Under  this  condition,  VH  ■  VS  5  0.  The  gas  dynamic  equations 

are : 

d&rrpAJ  +  l  dirrfiA.)  _ 

dl  ,  (,a) 

9 _  3((4,  +  4,)r) 

dtp  Qj  '  “  0  (lb) 

*  “  j1  -  2~  W  ~  ^)}“  (  U) 

where  A  ■*«'/««.  4,  —  mj «»,  A,  —  *»,/ a„  Am  -»  «■*  «/•„  —  p/p,,  p  —  p/p>, 

•  “«/«»»  —  t,TJm  —  ///m,  m—  \/(m—l),  w,  p,  p,  T , 

and  a  represent  relative  flow  velocity,  density,  pressure,  tem¬ 
perature  and  speed  of  sound,  respectively;  the  subscript  0  repre¬ 
sents  the  parameters  of  point  0  in  Figure  2(b);  $p  and  k  are  the 
isobaric  specific  beat  and  thermal  insulation  coefficient,  res¬ 
pectively.  All  the  linear  dimensions  used  are  the  ratios  with 
respect  to  the  width  B  in  the  l  direction  of  the  cascade;  t  is 
the  thickness  of  the  flowing  layer. 

From  Equations  (la)  and  (lb),  one  can  introduce  a  flow  func¬ 
tion  and  "position"  function  ♦  as  follows: 

Q^t/dtp  trfiAt,  df/tii  *  -  rpA,  (Id) 

Q<p/dt  m  4,,  df/dtp  —  (4,  +■  Am")r  (Je) 
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Figure  1.  Cascade  on  an 
arbitrary  stream  sheet 
of  revolution.  Refer¬ 
ence  lines  lj  and  V?  can 
be  selected  arbitrarily 
as  indicated  by  the 
dotted  lines. 


Please  note  that  the  presence  of  the  function  <f>  does  not  necess¬ 
arily  mean  that  the  flow  is  rotationless.  It  only  indicates 
that  the  component  of  the  absolute  degree  of  revolution  V  x  c 
perpendicular  to  the  flow  direction  is  zero  which  means  that  the 
flow  surface  coincides  with  the  revolution  surface. 


Figure  2(a)  Non-homoentropic  flow 
(the  state  point  N  can  move  in  the 
horizontal  direction  too)  K 


Figure  2(b)  Homoentropic 
flow 

1 — relative  stagnation  state 


2 .  Non-homoentropic  flow  (Figure  2(a) 

Under  the  condition  VH  =  0,  but  VS  K  0  and  D^/Dt  M  Q  (usually 
DS 

/Dt  is  used  to  approximately  consider  the  loss).  At  this  time, 
the  speed  of  sound  aQ  (and  T  )  at  point  0  in  Figure  2(a)  are  con¬ 
stant  in  the  entire  flow  field.  The  following  dimensionless  quan¬ 
tities  are  further  supplemented  in  the  following: 
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p*~p/pm*  fit  ““  p/ P*tt  S  “  S/e ,  *"  S/i 


The  subscript  00  indicates  that  the  parameters  are  for  point 
0Q  in  Figure  2(a).  Thus,  the  following  gas  dynamic  equations  [7] 
can  be  obtained:  a(r rfiA.r^  ,  8(rM.^^  _  0 

dl  dip 

dA,  _  dliA.  +  Am)Ll  _  =^L  rrp.%^ 

dq>  dl  K  difi 

fi  -  y  -  s»-  -  {1  _  [01*  -  A*)/(2m)]}-| 

A»  —  pexpC  —  AS) 

(p«/pl)e*p(  — AS)  —  1 

where  the  flow  function  ip  is  defined  by  the  following  equations: 

difi/dep  —  rrMi<<“4*'  “  f'fitA ,  1  (2d) 

difi/dl  —  -rfiA,*1-**  —  -tfitfij 


(2a) 

(2b) 

(2c) 


A  generalized  flow  function  0  can  be  defined  as: 

d9/dip  —  xrfiA,,  dB/dl  —  -rfi\  (2e) 

and  AS  *  mAS  =  S  -  Sqo  is  usually  given  in  the  form  of  A?  =  f(.l,4>). 

As  for  the  various  boundary  conditions,  they  will  be  given 
in  the  form  of  natural  boundary  conditions  in  the  variational 
principles  discussed  in  the  following  sections. 

II.  The  generalized  variational  principles  of  the  homoentroplc 
flow 


In  the  generalized  variational  principles,  the  variations  of 
various  gas  dynamic  parameters  are  independent.  From  this  the 
boundary  conditions  of  all  the  gas  dynamic  equations  of  cascades 
or  channels  on  an  arbitrary  stream  sheet  of  revolution  can  be 
automatically  derived  (including  the  continuity  conditions  on  the 
artificial  boundaries  E). 
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(A)  Generalized  variational  Principle  I 


The  necessary  condition  for  a  general  function  ^  which  is 
equivalent  to  the  solution  of  the  homoentropic  flow  of  a  cascade 
with  an  arbitrary  revolution  stream  sheet  is:  =  o. 


Here,  A„A9,<p,p  and  p  undergo  independent  variations.  If 
we  let  JK$.  -»«?/*)[  1- 1* < ?*->)]-  fi)  then> 

—  |||— *  -g|  -  ^  ^  gjr  ~  4  (.A*  “  —  n}trd<pdl 


+  f  (Am  +  Amir')(p  •  dS  +  (A  +  —  [</>(/)  Ip.}  •  dS 

J  J  W,) 

- ■;  (  ( A±  +  Aj  +  2/t.i,X«Ah  —  <!>y  —  •  dS 

+  4-  (  <*-  -  t+XA.  +  A.  +  2AJJ  •  dS  CM ) 

where  (A)  represents  the  area  of  the  region  where  solution  is 
sought  on  the  arbitrary  stream  sheet  of  revolution  Ci.e.,  the 
region  in  Figure  1);  (E)  represents  the  entire  arti¬ 

ficial  boundaries  (Figure  3);  Cc, }  represents  the  Inlet  surface 
j}*  _  -L 

and  the  outlet  surface  22*  (Figure  1);  (c2)  represents  the 
boundary  surfaces  13,  1 ' 3*  ,  ^  and  4*2*  in  the  periodic  region 

f 

in  the  front  and  back  of  the  cascade;  (c2)  represents  the  lower 
boundary  surface  T3  and  W;  (c^)  represents  the  S  and  p  surfaces 
along  the  profile  contours  of  the  blade;  (c)  represents  the  entire 
boundary,  i.e.,  (c^+^+c^);  the  subscript  +  and  -  represent  the 
the  right  and  left  side  of  the  integration  path,  respectively);  pr 
represents  the  predicted  quantity;  "  J:"  represents  the  upper 
boundary  surface  3*3'  and  4 '2 7 ;  "  T  "  represents  the  lower  bound¬ 
ary  surface  TJ  and  77.  In  addition,  A4'1>2  represent  the  differ¬ 
ences  of  <p  between  TJ  and  1 ' 3 '  and  T2  and  4  f2  f ,  respectively: 


I****  _ 

(rrpAi),jdip  (A< p  T*pr'C>9€nT 

*  sn^o p  itfck? 

£Translators  note:  pages  1-6 
X  *  top  and  *  bottom.  ] 
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Apparently,  is  just  equal  to  the  amount  of  gas  flow 

ejected  from  the  profile  contours  of  every  blade. 

Proof:  Taking  the  variation  of  (1-1)  and  using  the  Green 
theorem,  we  get: 

«*  -  H  {(aJ^«  - 1)  y  ■ ~  (7  ft  *■ ») 

+  4  CA,.„  —  AJWS  +  (4  -  <£*) 9A.JS 

» 

+  2*  Lp  +  ^  *>  </»>  -  A*.  4  «iiT»  •  dSf 

+  \  1<«  +  4^-X><4  —  x_)  +  (*„  -  ^„XM+  4-  »A.  )}  •  dS 

Key:  1,3, 5, 7 — upper;  2, 4, 6, 8 — lower 

Since  all  the  variations  are  completely  independent  and  arbitrary, 
therefore,  the  following  set  of  conditions  for  an  extremum  can 
be  obtained  of  6J^  =  0: 


Euler  equations: 

ajUg-fJUj  dA, 

dl  dtp  “ 

dip/dl  —  —  vfiA9t  d<f>/d<p  "•  rrpAi 
<*/«  4-  (A1  -  A’)/(2m)  -  1 
P-?'* 

a  Natural  boundary  conditions: 

<\  A^  —  [/!,(?)  1* 

Jdi C,  At,,  —  /It,,  and  <P±  “  <At  4-  AI1.1 


(extent  of  rotation) 
(equivalent  continuity 
equation) 

(energy  equation) 
(bomoentropic  relation) 


it  can  be  deduced  that  on  13  and  1'3'  Cas  well  as  on 
f?  and  4*2'),  all  the  gas  dynamic  parameters  are  equal 
(periodicity) . 

f&k  e,  <A,.,  —  [ <A,.,(0 ]„  (when  there  is  no  mass  transfer 

on  the  surface  of  the  blade,  then  ^  ,  =  given  constant) 

im  -  (/U),  and  «*,_  -  4>t  from  this  one  can  deduce 

that  all  the  gas  dynamic  parameters  are  continuous  on  Z. 


(I-2a) 


(I-2b) 
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It  has  already  assumed  that  q  *  pA  and  the  subscript  t 
indicates  the  tangential  component. 


It  is  obvious  here  that  from  this  generalized  variational 
principle  the  solution  to  the  homoentropic  linear  flow  of  a  cas¬ 
cade  satisfi^^e- all  the  boundary  condit^pns  and  is  conti¬ 

nuous  4tee  artificial  boundar^e^^^aan^b'/  efetainodi 

(B)  Generalized  variational  Principle  II 


The  necessary  condition  for  a  function  which  is  equiva¬ 
lent  to  the  homoentropic  flow  solution  of  a  cascade  with  an  arbi¬ 
trary  stream  sheet  of  revolution  to  have  an  extremum  is:  =  0 


Here  A^,  P  and  p  undergo  independent  variations 

Ju&Ai*  Ap,  4>,  p,  —  jj  -qJ~  +  y  ^  C.l*  +  A:.  +•  2AmAp) 

+  nl  rrdtpd!  —  (f ,)p ,4>rdS 

+  ~  I  ,  —  A£i.i)(qt  +  q*)  •  i ijrdS 

2  Jlfyi 

+  j  +  q-]  •  nriS  (Il-I) 


where  n,  t  represent  the  unit  vector  in  the  outside  normal  and 
tangential  directions  of  the  integration  path,  respectively; 

A<j>i>2  represent  the  differences  of  <J>  on  IT  and  1*3*  as  well  as  ^2 
and  4*2',  respectively: 

Jff&f  s 

^  ['■(A,  +  Aa)]iadq> 


Proof:  Taking  the  above  variations  and  from  6JII=0,  the  fol¬ 
lowing  set  of  necessary  conditions  for  an  extremum  can  be  obtained 


Euler  equations: 
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d(rr^i) 

di  +  d<j»  "" 
8<*>/di  —  Am  d*/d<J>  “  '<Af  +  A.) 
( 1\  +  CA»  -  A’)/C2«)  -  1 
K>'  *-0* 


(continuity  equation) 
(equivalent  rotational 
degree  function) 
(energy  equation) 
(homoentropic  relation) 


( II-2a) 


Natural  boundary  conditions: 

cifi  qi  —  U«Op)U 

—  q-r.,  and  **  “  +  &<*.. 

from  these  one  can  deduce  that  all  the  gas  dynamic 
parameters  are  equal  (periodicity) 

^OIK  <■_!  q»  —  (4.(0 !».  (where  there  is  no  mass  trans¬ 
fer  on  the  blade.,  qu=0) 

i  J£.  (*►>.-  W-j.  and  -  4*- 
from  the  above,  we  know  that  all  the  gas  dynamic 
parameters  are  continuous 


(II-2b ) 


It  is  apparent  that  from  the  generalized  variational  prin¬ 
ciple  II,  a  solution  to  homoentropic  flow  of  a  cascade  -which 
satisf-^s-  all  the  boundary  conditions  and  is  continuous 
artificial  boundar&s**aarr  be*" obtained. 


Ill .  Sub-generalized  variational  principle  family  of 
homoentropic  flow 

Prom  the  above  generalized  variational  principles,  a  new 
series  of  variational  principles  can  be  derived  (it  can  also  be 
called  the  sub-generalized  variational  principles).  The  principle 
used  here  is  based  on  that  if  the  extremum  conditions  of  a  var¬ 
iational  problem  are  applied  to  the  original  variational  problem 
as  constraints,  then  the  solution  of  the  extremum  would  not  be 
changed  [9].  The  following  examples  show: 
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(A)  The  sub-variational  principle  family  derived  from 

(a)  Substituting  rows  2  and  (4)  in  Equation  (-2a)  into  and 
■  eliminating  A1,  A«p  and  p?,  we  get  the  following  variational  principle: 


J'i  (*, 

J-  [(1 

2r‘j  1 V  r  dq>/ 

■♦(srn-t 

dt 

+ 

m  yfi  —  ^  Vrdopdl  ■+• 

+  A.  )4>dS 

+ 

.*-  A  *t) 

rfi  dn  / 

— 

i  L  <*• 

—  —  A4>,„ 

.)  {2AJ,  •  t  -  i.  [ 

(j 

£).♦< 

fit).]!- 

+ 

tL<*- 

-  *+>  {2AJ, 

t  l\ f  dp  A 

+ 

► 

(i  2). 

<i'-i ) 

1 

We  only  have  to  let  (assuming  ij;  and  p  undergo  variation 

independently)  to  obtain  the  first  and  third  rows  in  Equation  (I-2a) 
and  Equation  (I-2b). 


(b)  Substituting  the  third  and  fourth  rows  in  Equation  (I-2a) 
into  and  then  eliminating  p  and  $  we  get  another  subvariational 
principle : 


A„  A,) 


rr-n 


where  X=l-  (A2-^)/(2m),  represents  the  various  linear  integration 
terms  in  J-^ 


ft 

Since  J^=0  (let  A,  A^  vary  independently  with  respect  to  A<p) 
we  obtain  the  1st  and  2nd  rows  in  equation  (l-2a)  and  (l-2b). 


(c)  Using  the  2nd,  3rd  and  4th  rows  of  equation  (I-2a)  as  added 
constraints  to  J^,  we  obtain  the  following  variational  principle  III: 

-  jjj  F,  .  dq>dl  +  L%  (III-l ) 


where 


’-[* 


+  ♦  A.V  rx  •>  — 

« 


)  *  ♦  «±-» 


2m 


Y*  ♦  t.  I, 
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By  letting  *  0,  the  first  equation  of  (I-2a)  and 

Equation  (I-2b)  can  be  obtained. 

(B)  The  sub-generalized  variational  principle  family 
derived  from  Jjj 

Similar  to  the  variational  method  used  for  J^,  we  get: 


(a)  Sub-generalized  variational  principle 

y;,o.  h p) -|j  {f  |(g*y  +  (7  JJ)  -  ^ 

+  2  -*•-  [('*  l*).  *  (-•  I;).] 

+  2  [('  8*).  *  ('5  df).j  "S  <«  -'■> 


Prom  <$Jjj  *  0,  one  can  derive  the  first,  third  and  fourth 
rows  of  Equations  (II-2a)  and  (II-2b). 


(b)  Sub-generalized  variational  principle 

A„  A.)  -  f  (  p  +  4s 

\  A  •  ' 

+  m  \  X  —  1  —  1  (4J  +  {4,  +  \ti4tf.1'  +■  Lt>  Cll  1  ) 

where  represents  the  linear  integration  terms  in  J^; 

—  _  vm 
P  =  X 

Prom  <SJjj  *  0,  the  first  and  second  rows  of  (II-2a)  and 
(II-2b)  can  be  derived. 


Cc)  Variational  principle  IV. 


Using  the  second, 

a+/*CrrZnu  ' 


_  _ 

as  Uie-  constraints 


to  J 


third  and  fourth  rows  of  Equation  (JI-2a) 
XI,  we  get 


L2 


where 


y.v<*) 


(IV-O 


H 


G  id tpdi  -■  ^ ii 


C, 


2  +  (1 2*Y 

2 a.  a^ii 

l\S/  ) 

+  \r  d<p) 

r  d?JJ 

Let  <SJjy  ®  Oj  then  the  first  Equation  in  (JI-2a)  and  Equa¬ 
tion  Cll-2b)  can  be  obtained. 

Similarly,  more  sub-generalized  variational  principles  can 
be  derived  from  and  Jjj*  They  will  not  be  further  discussed 
here . 


IV.  Dual  extremum  principle  of  homoentroplc  flow 

Now  we  are  going  to  prove  that,  under  subsonic  conditions, 
the  variational  principles  III  and  IV  discussed  above  all  belong 
to  the  extremum  principle.  Furthermore,  they  form  a  pair  of  dual 
extremum  principles. 

CA)  Extremum  (minimum)  principle  V 


The  homoentroplc  subsonic  flow  solution  of  a  cascade  on  an 
arbitrary  stream  sheet  of  revolution  makes  the  general  function 

J 1 1 1 )  to— nova  a  minimum. 


Proof:  It  can  be  proven  [1]  that  when  6Jiy  _  Q  there  is 


*VniO) 


-11 


1&- 


/// 


TA(?  -  A*) 


d<?4l 


It  is  not  difficult  to  see  that  the  condition  to  to  be 
positive  is  5*  —  A* >o  .  Therefore,  we  know  that  at  subsonic 
speed,  the  function  has  a  minimum. 


r 


(B)  Extremum  (maximum)  principle  VI 

The  homoentropic  subsonic  flow  solution  of  a  cascade  on  an 
arbitrary  stream  sheet  of  revolution  makes  the  function  ( 4> ) 
to  have  a  maximum. 


The  proof  is  similar  to  the  above.  It  will  not  be  repeated. 
(C)  Dual  extremum  principle  VII 

a)  channels  on  an  arbitrary  stream  sheet  of  revolution. 


Under  this  condition  the  boundary  sections  13  and  3T  shown 

in  Figure  1  do  not  exist.  The  distribution  of  cf>  along  c^  can  be 

directly  given.  Therefore,  the  term  (  is  no  longer  necessary 

Ju,* 

in  JjjjCiIO.  The  solutions  corresponding  to  the  extremum  Cexpressed 
by  using  the  subscript  ex)  can  be  used  to  write  Equations  (III-l) 
and  (IV-1).  They  are 


+  (A*  +  A»A,)pj  xrdcpdl 


(«) 


Applying  Gauss  theory  and  also  considering  Equations  Cla) 
and  tie),  we  can  write 

•  nrdS  — 

“  (j  («1  *  rdydl  «  ((  {p(A*  -4-  A«A,))„r  rd  pdl 

(A) 

Substituting  the  above  into  Equation  (.8)  and  then  comparing  with 
Equation  (a),  we  obtain  the  following  important  relation: 

7lll(<Au)  “  J IvC “  /„ 


This  is  the  dual  extremum  principle  for  channels  on  an 


arbitrary  stream  sheet  of  revolution:  Not  only  that,  the  extre¬ 
mum  solutions  4>ex  and  4>gx  make  the  general  functions  Jjjj  and 
have  a  maximum  and  a  maximum,  respectively,  but  also  the  two 
extrema  are  equal. 

Generally  speaking,  JI];IOP)  ) . 

b)  cascade  on  an  arbitrary  stream  sheet  of  revolution. 

Using  treatment  similar  to  the  one  described  above  and 
selecting  the  stream  lines  at  stationary  points  in  the  front  and 
rear  of  the  blade  contours  as  the  boundaries  IT  and  FT ,  we  get 


where 


ZlllC*..)  “  Jw(<pn)  +  /« 


CVII-O 


This  is  then  the  dual  exrremum  principle  for  a  cascade  on 
an  arbitrary  stream  sheet  of  revolution.  Not  only  that,  the 
extremum  solution  ipQx  and  <J>ex  make  the  general  functions  J-j-jj  and 
Jjy  to  have  a  minimum  and  a  maximum  respectively,  but  also  there 
exist  a  relation  between  the  two  extrema  as  described  by  Equation 
CVII-1). 

V.  Generalized  variational  principles  and  variational  principles 
of  non-homoentroplc  flow 

Here  the  distribution  of  entropy  S  *  is  given  (.Figure 

2  (.a ) ) . 

(A )  Generalized  variational  prlnc iple  (VIII) 


The  condition  under  which  the  general  function  Jy  which  is 
equivalent  to  the  nonhomoentropic  flow  solution  of  a  cascade  on 
an  arbitrary  stream  sheet  of  revolution  is:  6Jy  =  0. 
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Here  ,  A„  A9t  <fi,  p,  and  pQ  all  undergo  variation  independently. 
MA„  A)  -  jj  {£  “  (-^yA')|f  ~  ft*  ~  Al )  -  i?}  Trdq>dl 


where 


+  +  —  i<koi»«}  •  a 

—  yj^  (4j.  +  4y  +  2/1.4)(^  —  4?  —  &4>iu)  • 

+  y  Jtx,  <*-  -  *»X4+  +  4.  +  2A.i^  •  dS  ( v-i  > 

0(*,  A,  ;)  -  "•  {  j»  [l  +  AJ  -  in  CM- o]  -  A  } 


Proof:  The  proof  is  similar  to  the  one  used  in  the  treatment 
of  generalized  variational  principle  I.  Prom  6JV  »  0,  we  can 
obtain  the  entire  set  of  Euler  equations  and  natural  boundary 
conditions.  With  the  exception  that  the  degree  of  revolution 
equation  and  the  homoentropic  relation  are  replaced  by  (2b)  and 
(2c) respectively,  and  p  and  p  are  replaced  by  pQ  and  pQ, 
respectively,  the  rest  are  identical  to  Equations  (I-2a)  and  (I-2b). 
It  can  be  seen  that  from  the  generalized  variational  principle  VIII 
we  can  derive  a  nonhomoentropic  flow  solution  of  cascade  which 
satisfies  all  the  boundary  conditions  and  is  continuous  on  the 
artificial  boundaries. 


(B)  Sub-generalized  variational  principle  family 

Using  the  method  similar  to  the  one  used  in  Section  III,  we 
can  derive  another  family  of  sub-generalized  variational  principle. 
Only  one  of  them  is  listed  here. 


Variational  principle  IX 

hiU>)  -  jj  '  i'  *  A‘  *  KA'  ^ 

X  [l  -  CA*  -  A*)]“  d<pdl  +  LiU 


(Vt-1) 
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l 


f 


where  Ljjj  represents  all  the  linear  integration  terms  in  Jy. 
Prom  SJyi  *  0,  we  can  obtain  the  degree  of  revolution  Equation 
(2b)  and  the  entire  (natural)  boundary  conditions. 


(C)  The  variational  principle  X  expressed  by  the  general¬ 
ized  flow  function  &  ,  The  condition  for  the  following  function 
Jyn  which  is  equivalent  to  the  nonhomoentropic  flow  solution  of 
a  cascade  on  an  arbitrary  stream  sheet  of  revolution  is  6JVII  *  0. 


.  i 

/vu(S)  “*  (]  ^ie*P (—•&£)d<pdl  +  —  AS)(A„  +  A.i,)8  •  dS 

+  fu)cxp (-ASXA  +  -  (SCOlp.)  •  dS 

-  —  (  ,  exp(-A5)M±  +  A?  +  IAJJQBm  -  9y  -  A©,.,)  •  dS 

+  —  (  «p(-ASX«-  -  0+X/U  +  A-  +  ■  dS 

2  .)<*> 


Cvii— i ) 


From  *  0,  we  can  get  Equation  (2b)  and  all  Cnatural) 

boundary  conditions 


VI .  Regarding  the  actual  application  of  the  variational  principle 


The  key  here  is  the  treatment  of  the  various  complicated  bound¬ 
ary  conditions  f  the  cascade  problem.  We  primarily  rely  on  taking 
full  advantage  of  the  unique  functions  of  the  natural  boundary  con¬ 
ditions  and  the  artificial  boundaries.  Now  we  are  going  to  brief¬ 
ly  explain  the  situations  separately. 


1^In  JVII>  let  us  assume  that  S  *  S(_l,<|0  is  given.  In  practice, 
it  ismore  convenient  to  make  s-s(i.n)  .  But  we  agree  that 
it  does  not  undergo  variation  which  means  that  6 CAS)  *  0.  Under 
such  conditions,  the  variational  principle  X  is  still  valid. 
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In  all  the  above  variational  principles,  all  the  boundary 
conditions^  have  been  transformed  into  natural  boundary  conditions. 
Therefore,  they’  do  not  have  to  be  considered  in  seeking  for  the 
solution.  Despite  this  fact,  in  the  direct  variational  methods 
(such  as  Ritz  method  and  Kantorvich  method,  etc.),  it  is  still 
more  advantagenous  if  the  coordinate  function  can  satisfy  some 
of  the  natural  boundary  conditions.  This  would  improve  the 
accuracy  of  the  solution.  However,  it  is  not  practically  poss¬ 
ible  to  choose  a  coordinate  function  beforehand  which  will  satisfy 
all  the  boundary  conditions  in  the  region  wtx-n  a  solution  is 
being  sought.  It  appears  that  the  most  effective  and  simple 
method  is  to  use  the  artificial  boundar-^with  high  flexibility  to 
choose  the  coordinate  function  by  the  region.  That  is  to  analyze 
the  characteristics  of  difference  in  e?.ci  region  then  first  to 
separate  them  and  followed  by  solving  them.  For  example,  in 
Figure  1,  the  characteristics  of  the  boundary  conditions  between 
zones  I,  inJ^^Ere^* drastically  different.  Therefore^it 
more  proper  to  c'jfoose  the  straight  line  33Y  and  44  '/as  t^e^rti- 
ficial  boundary  Z.  Henceforth,  in  regions  I  and  III,  the  coordi¬ 
nate  function  only  must  satisfy  the  periodicity  requirement  and  in 
region  II,  it  only  must  satisfy  the  blade  surface  conditions.  In 
the  meantime,  because  the  linear  integration  terms  on  the  Z  surface 
have  already  been  added  to  all  the  variational  principles,  the 
approximate  solution  combined  by  all  the  regions  would  then  auto¬ 
matically  (although  it  is  an  approximation,  yet  it  is  the  highest 
limit)  satisfy  the  continuity  condition  on  Z. 

The  artificial  boundaries  can  be  used  to  separately  process 
the  various  boundary  sections  with  different  geometric  character¬ 
istics.  For  example,  in  Figure  3  aa '  and  64 T  divided  the  surface 
of  the  blade  into  sections  2a,  ab,  PT,  3 'a'  and  a ' H '  with  differ¬ 
ent  geometric  properties  (arcs  of  various  radius,  sections  of  a 
parabola,  etc.)  and  further  divided  region  II  into  sub-regions 
II,  IIj  and  II^.  Therefore,  the  blade  surface  curve  in  each  sub- 


18 


region  can  easily  be  expressed  or  approximated  by  simple  analy¬ 
tical  equations.  It  is  then  not  difficult  to  structure  the  coor¬ 
dinate  function  which  satisfies  the  boundary  conditions  of  the 
blade  surface.  Apparently,  for  the  simplification  in  the  calcula 
tion  of  various  integrations  j  ,  the  artificial  boundaries 
should  be  straight  lines. 


Figure  3-  The  method  to  select 

the  artificial  boundaries 

Key:  i)  separate  the  regions  where 
the  characteristics  of  the 
boundary  conditions  are  diff¬ 
erent  (such  as  33'  and  W1"). 
ii)  separate  the  boundary  sec¬ 
tions  with  different  geomet¬ 
ric  characteristics  (.such  as 
aa  ’  and  ** 


In  the  finite  element  method,  the  boundary  of  each  element 
can  be  used  as  the  artificial  boundary  E.  Therefore,  if  more 
undetermined  parameters  are  added  to  the  implantation  function  of 
each  element,  through  the  action  of  the  I  terms  in  all  the 

wit 

general  functions,  it  will  automatically  attain  an  optimal  con¬ 
tinuity  (i.e.,  miscibility)  on  the  boundary  E  of  every  element. 

The  convergence  and  accuracy  of  the  solution  can  then  be  assured. 

In  either  the  finite  element  method  or  the  direct  variational 
method,  after  an  approximate  function  is  substituted  into  the 
general  function,  the  solution  can  be  obtained  directly  using 
various  advanced  mathematical  calculation  methods. 

VII .  Conclusions 


1.  In  the  above  three  families  of  variational  principles 
and  dual  extremum  principles  of  homoentropic  and  nonhomoentropic 
flow  for  the  direct  problem  of  cascades  and  channels  on  an 
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arbitrary  stream  sheet  of  revolution  were  derived.  The  major 
characteristics  are:  a)  it  considered  the  conditions  under 
which  gas  injection  or  suction  along  the  surface  of  the  blade 
exists  in  order  to  satisfy  the  requirements  of  the  newly  developed 
area  such  as  boundary  layer  control,  air  film  type  or  prespira- 
tion  type  of  cooling,  etc.;  b)  it  takes  full  advantage  of  the 
special  functions  of  the  natural  boundary  conditions  and  the  arti¬ 
ficial  boundaries  to  simplify  the  treatment  of  the  boundary  con¬ 
ditions.  Thus,  the  practical  value  of  these  variational  principles 
can  be  raised. 

2.  In  addition  to  the  fact  that  each  generalized  variation¬ 
al  principle  can  be  used  to  directly  solve  problems  (through  the 
finite  element  method,  variational-difference  method  [10]  or  the 
direct  variational  method.  The  major  advantage  is  that  the  dis¬ 
tribution  of  flow  velocity  does  not  have  to  be  derived  from  partial 
differentiation  with  respect  to  (or  <f>,  or  .91  so  that  better 
accuracy  can  be  expected.  The  disadvantage  is  that  the  solution 

is  not  uniquely  determined.  Therefore,  it  is  more  difficult  to 
analyze  convergence  and  error),  another  important  feature  is  that 
a  series  of  new  sub-generalized  variational  principle  families  can 
be  derived  from  them  and  thus  providing  an  wide  degree  of  freedom 
to  choose  from  in  its  application. 

3.  Using  the  variational  principles  in  this  paper  which  con¬ 
sidered  the  mass  transfer  on  the  surface  of  the  blade,  the  condi¬ 
tion  under  which  suction  on  the  blade  surface  exists  can  be 
treated  rigorously.  For  the  condition  of  gas  injection  on  the  sur¬ 
face  of  the  blade,  it  was  assumed  that  the  H  and  k,  etc.,  of  the 
injected  gas  (except  §)  are  the  same  as  the  main  gas  flow.  If  the 
H  is  different  for  the  injected  gas,  a  corresponding  variational 
principle  can  also  be  derived.  This  will  be  discussed  elsewhere. 
This  method  can  also  be  used  to  simulate  the  repulsive  reaction 

of  the  boundary  layer  of  the  blade  surface  and  tail  trace. 
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4.  The  various  variational  principles  introduced  in  this 
paper  are  only  applicable  to  subsonic  flows  and  non-turbulent 
transonic  flows  (the  partial  supersonic  section  included  must 
undergo  special  treatment,  refer  to  reference  [11]). 

5.  It  can  be  proven  that  the  physical  meaning  of  the  surface 

integration  part  of  Jjj-j-  aBd  as  well  as  other  general  func¬ 
tions  J^,  J^,  J" ,  and  Jv  derived  from  the  same  family  is  the  diff¬ 
erence  of  kinetic  energy  K{—cx/l')  and  potential  energy  (internal 
energy)  n  (=CvT)  of  the  flowing  matter  in  the  region  where  the 
solution  is  obtained.  Considering  further  the  generalized  enthalpy 
H:  —  n  +  H)f*iV  it  can  be  seen  that  this  family  of  the 

(4) 

variational  principles  is  very  closely  related  to  the  Hamiltonian 
principle.  Furthermore,  the  physical  meaning  of  the  surface  inte¬ 
gral  part  of  Jpy  and  the  other  general  x’unctions  of  the  same 
family,  such  as  Jjjj  Jj-j  and  is  the  kinetic  power  of  the  total 
flowing  material  in  the  region  where  the  solution  is  obtained. 

6.  It  goes  without  saying  that  the  entire  variational  prin¬ 
ciples  discussed  In  this  paper  can  be  applied  to  treat  the  special 
case  of  planar  flow  conditions. 

7.  This  paper  considered  the  effect  of  viscosity  loss  in 

approximation  by  properly  estimating  the  distribution  of  entropy 
SCI, 'll).  If  a  more  perfect  viscosity  model  is  adopted,  which  is  to 
add  a  a  viscosity  term  to  the  right  of  Equation  (2b)  and 

then  to  consider  it  as  a  known  quantity  (it  can  be  obtained  by  an 
iteration  method),  it  only  requires  a  slight  modification  on  the 
nonhomoentropic  variational  princioles  discussed  In  the  paper  (by 
adding  a  surface  integral  term  •  Here  —  }+/ 

is  the  circumferential  component  of  the  dimensionless  viscosity. 
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CALCULATION  OF  UNSTARTED-CHOKED  REGIME  OF  SUPERSONIC  AND  TRAN¬ 
SONIC  COMPRESSOR  CASCADES* 

Cheng  Jia-gang 

(Institute  of  Mechanics,  Academia  Sinica) 

ABSTRACT 


This  paper  studiei  unstarted-choked  regime  of  supersonic  and  transonic  cascade, 
such  regime  is  important  for  the  starting  process  and  offdeaign  regime. 

A  semiempirical  method  is  given  here,  which  allows  for  calculating  operating 
incidence  of  this  regime.  Under  the  condition  of  downstream  throttle,  or  at  various 
the  axial  velocity  density  ratio,  this  method  also  gives  useful  results.  Other  perfor¬ 
mance  parameters,  such  as  total  pressure  loss,  pressure  ratio  and  flow  deviation  angle 
may  be  obtained  simubtaneously. 

In  comparison  with  other  works1’"*1,  the  method  presented  here  considers  more 
variation  of  condition.  The  calculating  results  obtained  -show  this  method  is  reliable. 
Because  this  method  is  more  simple  and  convenience,  it  may  be  adopted  to  analyse 
the  aerodynamic  performance  of  cascades  and  the  influence  of  geometric  parameters 
and  dynamic  parameters  on  cascade  performance. 


This  paper  studied  the  unstarted-choked  regime  of  supersonic 
and  transonic  cascades.  Under  conditions  of  downstream  throttle 
and  varying  axial  velocity  density  ratio,  useful  results  could 
also  be  obtained.  In  addition,  other  corresponding  performance 

v  AV  -■  .  .•  f- 

parameters  of  the  cascade,  such  as  total  pressure  loss4  pressure 
ratio  and  flow  deviation  angle  can  also  be  obtained.  In  compar¬ 
ison  with  other  works  reported  by  foreign  researchers ,  this  method 
may  be  used  to  take  more  variations  of  condition  into  considera¬ 
tion.  The  results  of  our  calculation,  after  the  comparison,  indi¬ 
cated  that  this  method  Is  reliable.  Since  this  method  is  rela¬ 
tively  simple,  it  may  be  used  to  analyze  the  performance  charac¬ 
teristics  of  the  cascade.  Furthermore,  it  may  be  used  to  analyze 
the  performance  characteristics  of  the  cascade  under  conditions 
that  geometric  and  dynamic  parameters  vary. 

1 - 

Received  on  April  3,  1978. 
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SYMBOLS 


M  Mach  number 

K  experimental  constant 

w  gas  flow  velocity 

a  speed  of  sound 

0  gas  flow  angle  (.in  reference 

to  the  axis  of  the  blade) 

b  thickness  of  the  flow  layer 

p(p°)  static  pressure  of  the 
gas  flow  (total  pressure) 

p  density  of  the  gas  flow 

L  hydrodynamic  factor  of  the 

cascade 

V  Reynolds  factor  of  the 

cascade 

U  Mach  factor  of  the  cascade 

a  density  of  the  cascade 

8  installation  angle 

s 

<i>  deviation  angle 

R  arc  length  Reynolds 

number 


J  efficiency  correction  factor 

ip  thickness  at  the  rear  end  of 
the  blade 

Q  axial  ¥8 1  aciifc y  density  ratio 

S^/t  compression  ratio  at  the 
throat 

A2*/t  downstream  throttle  ratio 

Y  specific  heat  ratio 

q(M)  critical  compression  ratio 

n  subsonic  cascade  varying 
efficiency 

c  arc  length 

t  nodal  distance 

oT  coefficient  of  total  pressure 

loss 

Pp/p,  cascade  pressure  increase 
xratio 

A0*81-82  gas  flow  deviation  angle 

T  temperature  of  the  gas  flow 

Q/Qmov.  flux  ratio 
max 

Subscripts : 

1.  Flow  far  in  front  of  the 
cascade 

2.  Flow  way  downstream  from  the 
cascade 

No  subscript: 

Flow  behind  the  shock  wave  in  the 
duct 
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There  are  three  types  of  inlet  flow  patterns  for  supersonic 
cascades:  unstarted  flow  pattern  Cor  mixed  type),  started  flow 
pattern  (supersonic  type)  and  saturated  flow  type  (.supersonic 
axial  velocity  type).  The  characteristic  of  the  first  type  is 
that  there  is  a  quasi-normal  shock  which  is  detached  from  the 
object  in  the  inlet  region.  The  second  flow  pattern  is  the  appear¬ 
ance  of  an  oblique  shock  wave  which  is  attached  to  the  object  in 
the  inlet  region.  It  is  supersonic  In  the  duct  and  in  a  started 
condition.  The  back  pressure  from  downstream  cannot  influence  the 
flow  pattern  of  the  inlet  area.  '  The  saturated  flow  pattern  appears 
when  the  Mach  number  of  the  incoming  flow  is  even  higher  (.for 
example  M^  >  2).  At  this  time,  the  velocity  component  along  the 
axis  of  the  cascade  becomes  supersonic.  The  shock  wave  and  Mach 
wave  are  limited  to  the  inside  of  the  duct  of  the  cascade.  During 
the  starting  process  of  the  cascade,  when  varies  from  lower  to 
higher  number,  these  three  types  of  flow  pattern  can  appear  In 
sequence.  The  superiority  of  the  saturated  flow  pattern  is  only 
limited  to  the  supersonic  aerospace  vehicles  which  has  not  been 
studied  extensively  at  this  moment. 

In  a  compressor  with  shock  wave  inside  the  rotor,  as  soon  as 
after  the  entire  blade  height  reaches  supersonic  condition,  the 
characteristic  will  appear  to  be  very  steep  as  shown  in  Figure  1. 
The  vertical  section  represents  the  unstarted  choked  flow  pattern. 
Its  upper  limit  of  Mach  number  (=1.2-1  5)  is  determined  by  the  area 
of  the  throat  region  and  the  dullness  of  the  front  fringe.  The 
curved  section  represents  the  unstarted  unchoked  flow  pattern.  Its 
study  will  assist  In  the  prediction  of  loss  of  velocity.  After  M-^ 
has  exceeded  the  upper  limit,  the  cascade  then  moves  into  the 
started  working  condition.  The  characteristic  steepness  is  created 
by  the  sole  attack  angle  Cor  the  periodicity  of  the  flow  limits 
the  flux).  Apparent  in  the  characteristics  of  variation  of  working 
conditions  in  a  compressor,  there  are  many  situations  Cespecially 
under  low  revolution  speed  and  high  back  pressure  conditions)  which 
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belong  to  the  unstarted  choked  flow  pattern. 


Figure  1.  The  characteristics 
of  a  supersonic  compressor 

Key:  1 — curved  section 

2 — vertical  section 


Presently,  more  working  conditions  are  designed  to  use  the 
started  flow  pattern.  However,  the  unstarted  choked  flow  pattern 
has  also  been  used  to  design  the  working  condition  of  the  first 
order  stator  (.to  prevent  the  propagation  of  noise  forward).  In 
addition,  it  was  proved  in  [1]  that  the  use  of  the  unstarted  in¬ 
tense  oblique  shock  wave  series  design  for  rotors  is  more  superior 
than  using  the  multiple  arc  cascade  at  high  pressure  ratios.  This 
paper  studied  the  unstarted  choked  working  condition  of  supersonic 
cascades.  In  the  judgment  of  the  operating  attack  angle  of  this 
type  of  flow  pattern,  there  have  been  many  researchers  abroad 
studying  this  problem  using  semi-empirical  methods  [2-5].  This 
paper  provided  a  relatively  simple  method  to  calculate  the  operat¬ 
ing  attack  angle  of  the  choked  flow  pattern  encountered  during  the 
starting  process  or  the  variation  of  working  conditions.  It  also 
considered  the  effect  of  the  back  pressure  from  downstream  and  the 
axial  density  ratio  on  the  operating  attack  angle. 

In  order  to  calculate  the  attack  angle  8^,  it  is  necessary 
to  study  other  characteristics  of  this  flow  pattern,  such  as  loss, 
pressure  ratio,  etc.  Owing  to  the  complexity  of  the  problem  very 
often  people  would  adopt  various  postulated  flow  model  and/or  semi- 
empirical  calculation  methods  which  Ignore  the  position  and  shape 
of  the  shock  wave  as  well  as  its  connection  with  the  characteristics 
of  the  boundary  layer  in  order  to  study  the  main  factors  Influencing 
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the  loss.  They  have  provided  very  valuable  information  [6,7]  in 
the  areas  of  research,  design  and  experimentation.  This  paper 
referred  to  the  treatment  used  in  [73  to  offer  a  simple  model  of 
cascade  loss  to  calculate  other  properties  of  the  cascade. 

I .  ANALYSIS 

Because  the  characteristic  curve  of  supersonic  and  transonic 
compressors  is  steep,  it  is  approximately  an  equi-mass  flow  rate 
line,  the  correct  estimation  of  its  operating  attack  angle  B-l  is 
very  important.  A  1°  difference  in  B-^  can  result  to  a  5%  differ¬ 
ence  in  flux.  One  of  the  reasons  why  the  early  stage  failure  in 
the  development  of  the  supersonic  compressor  is  the  lack  of 
accuracy  in  calculating  the  flow  angle.  It  was  impossible  to 
design  the  blade  profile  with  small  loss.  Thus,  the  efficiency 
was  very  low. 

The  flow  of  a  planar  cascade  is  basically  two-dimensional. 
When  M1  of  the  incoming  flow  is  higher,  shock  ±s  created.  The 
starting  shock  wave  originates  as  a  supersonic  bag  on  the  back  of 
the  blade  and  begins  to  develop  into  a  quasi-normal  shock  wave 
across  the  entire  duct  with  increasing  M^.  Based  on  the  suggestion 
by  Perri  [2],  the  intensity  of  the  shock  wave  in  the  duct  is  higher 
than  the  heat  wave.  Therefore,  in  the  one-dimensional  analysis, 
it  only  considers  the  proper  making  of  assumptions  of  intensity 
and  position  with  respect  to  the  shock  wave  in  the  duct  (it  is 
also  called  the  starting  shock  wave  in  the  starting  process).  It 
does  not  care  about  the  actual  shape  of  the  shock  wave.  The  heat 
wave  is  neglected  completely.  Many  authors  assumed  that  the  un¬ 
started  flow  pattern  is  the  normal  shock  wave  inlet  model.  There¬ 
fore,  the  gas  flow  is  not  deflected  after  the  shock  wave  and  the 
intensity  of  the  wave  is  only  affected  by  the  incoming  flow 
which  is  realtively  crude.  Actually,  In  addition  to  the  effect  of 
M-p  the  shock  wave  In  the  duct  is  also  Influenced  by  many  factors. 


such  as  the  reverse  pressure.  Furthermore,  the  flow  behind  the 
wave  Is  deflected. 

For  simplicity,  let  us  assume  the  flow  is  one-dimensional 
and  the  shock  wave  in  the  duct  deflects  the  incoming  flow  from  6-^ 
direction  to  3  direction.  The  proper  selection  of  the  6  value  can 
improve  the  accuracy  of  the  analysis.  Apparently,  the  choice  of  8 
is  affecting  the  outcome  of  the  calculation  significantly.  If  we 
choose  8  =  is  then  equivalent  to  the  normal  shock  wave  case. 

Chauvin  [4]  suggested  that  the  average  direction  of  the  duct  at 
the  inlet  of  the  cascade  be  8.  However,  Fabri  [3]  recommended  the 
direction  of  zero  rising  force  of  the  cascade  to  be  8.  In  this 
paper,  however,  we  used  the  geometric  average  direction  of  the  gas 
flow  directions  at  the  inlet  and  the  outlet  as  the  direction  of 
gas  flow  after  the  shock  wave  in  the  duct. 

For  compressible  viscous  gases,  average  geometric  direction 
and  zero  rising  force  direction  are  different.  Using  the  former 
to  substitute  the  latter  makes  it  much  simpler  in  the  treatment. 

Only  that  the  gas  flow  direction  at  the  outlet  of  the  cascade  62 
must  still  be  known.  In  the  position-flow  calculation  of  a  finite 
thickness  rounded  rear  fringe  cascade,  02  and  the  ring  quantity  are 
In  principle  uncertain.  This  is  also  the  common  drawback  of  all 
position  flow  calculation  methods.  Theoretically,  we  can  use  the 
rear  pressure  region  to  replace  the  rear  stationary  print  to  per¬ 
form  a  calculation  to  compensate  for  viscosity  to  solve  the  problem. 
However,  it  is  still  hard  to  solve  the  flow  in  separated  regions. 
Therefore,  many  research  work  switched  to  the  use  of  comprehensive 
data  obtained  from  cascade  experiments  to  seek  for  solution.  They 
performed  statistical  analysis  on  the  experimental  results  obtained 
from  various  types  of  wind  tunnel  and  various  subsonic  cascade 
wind  tunnels  with  different  geometric  parameters  and  dynamic  para- 

r —  CO^ 

meters  and  found  that  r  and  cos  2  has  a  dependence  relationship 
close  to  that  of  a  hyperbolic  curve  [8]  (Figure  2).  In  the  same 
figure  some  experimental  data  of  supersonic  cascade  and  supersonic 
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rotor  are  also  plotted.  Therefore,  in  this  paper  we  assume  that 
Q  cos  81/cos  82  *  K  is  the  supplementary  equation  for  the  gas  flow 
angle  at  the  exit.  The  value  of  K  should  vary  depending  on  the 
geometric  and  dynamic  parameters  of  the  cascade.  It  is  most  appro¬ 
priate  either  to  determine  it  experimentally  or  to  choose  its  value 
based  on  lots  of  statistical  data.  As  for  the  data  gathered  in 
Figure  2,  it  is  more  suitable  to  let  K  =  0.75-0.8.  The  physical 
meaning  of  K  can  be  considered  as  the  flow  density  ratio  at  the 
inlet  and  the  outlet.  It  is  apparent  that  in  a  compressor  cascade 
the  variation  of  pg/p-^  and  w2/w1  has  a  mutually  weakening  affect. 
Their  products  are  not  very  sensitive  to  the  variations  of  other 
parameters . 


Figure  2.  Relation  between  measured  flow 
angle  of  ccrpressor  cascade  and  Q. 


— 

mi 

NACA06J 

- ' 

u,<0.t 

— * 

t»i 

Cl 

9  —  5 

U,<9 1 

— — 

m 

Cl 

9-2 

M,«l 

m 

8TB 

9-3 

U,  -  1.56 

• 

poj 

MCA 

9St2 

W,  -  1.5 

* 

Ml 

Bator 

— 

U,  -  1.2—1.54 

In  C9],  some  crude  non-compressible  analysis  was  given  to 
this  hyperbolic  relationship  obtained  based  on  experimental  statis¬ 
tics.  That  paper  assumed  that  separation  occurred  at  a  certain 
point  on  the  back  of  the  blade.  The  local  diffusion  factor  is 
D„_  =  W...,  /Vf  Because  the  separation  region  remains  to  be 

sep  sep  max 

isobaric,  that  is  Dsep  is  invariant,  the  velocity  on  the  outside 
of  the  separation  region  is  W2  =  Wsep.  It  was  further  assumed  that 

the  main  flow  which  did  not  suffer  velocity  loss  could  still  be 

* 

deflected  to  the  nominal  02  angle  direction,  then  we  get: 


For  a  given  cascade,  is  approximately  invariant  and 

Dsep  *s  a  constant*  Therefore,  the  right  hand  side  of  the  above 
equation  Is  approximately  a  constant. 

Because  the  pressure  Increase  of  a  supersonic  and  transonic 
cascade  is  relatively  larger,  the  effect  of  compressibility  on  the 
average  geometric  direction  should  be  considered.  The  subscripts  u 
and  a  are  used  to  represent  the  facial  and  axial  components  of  the 
gas  flow  velocity  in  the  cascade,  respectively.  Then  after  dividing 

*"  W,  1  +  («>*/ ]/2 

wujwu  —  Qpjp„  ve  get 

when  0—1,  and  p^/ p2  .=  this  equation  is  reduced  to  tg|5f  =  CtgB^'*’ 

+  tgB2)/2.  It  is  obvious  that  in  order  to  solve  for  B-^,  we  must 
consider  other  characteristics  of  the  cascade,  such  as  the  cascade 
density,  etc.  To  calculate  the  cascade  loss,  it  is  not  sufficient 
to  calculate  the  shock  wave  in  the  duct  term  alone.  The  model 
adopted  by  this  paper,  in  general,  was  that  the  unstarted  supersonic 
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*shock  wave  in  the  duct  deflected  the  incoming  flow  to  the  aver¬ 
age  geometric  direction  of  the  cascade  and  stowed  it  down  to  sub¬ 
sonic.  The  flow  behind  the  wave  traveled  through  the  subsequent 
subsonic  compressors  in  a  one-dimensional  flow  manner.  The  opti¬ 
mal  variational  efficiency  of  compression  nQ  is  related  to  the 
hydrodynamic  factor  of  the  cascade  L,  Reynolds  factor' V  and  Mach 
number  factor  U  based  on  the  dimensional  analysis  given  in  £10], 
Other  losses  not  accounted  for,  such  as  shock  wave  boundary  layer 
interference,  tail  trace  mixing,  etc.,  are  represented  by  the  effi¬ 
ciency  correction  factor  f  obtained  based  on  the  statistical  treat¬ 
ment  of  experimental  data.  In  the  situation  that  experimental  data 
are  lacking,  the  results  listed  in  [7J  will  be  used  temporarily. 

Through  the  above  analysis,  we  can  briefly  describe  the  assump 
tions  made  in  our  calculation  as  follows: 

i.  The  flow  passing  through  the  cascade  is  basically  one¬ 
dimensional.  With  no  flow  restriction,  the  throat  is  at  the  same 
state . 

ii.  Neglecting  the  heat  wave,  we  assume  that  the  shock  wave 
in  the  duct  deflects  the  incoming  flow  to  the  average  geometric 
direction  of  the  gas  flow  at  the  inlet  and  outlet. 

iii.  There  exists  the  experimental  statistical  rule 
Qco*fijco»fi,  —  K  and  K  Is  determined  experimentally. 

It  The  subsonic  cascade  loss  behind  the  shock  wave  in  the  duct 
is  mainly  due  to  the  compression  loss  nQ  and  separation  loss  f. 

II.  DERIVATION  OF  EQUATIONS 

Let  M,  6,  p,  p°  and  b  represent  the  Mach  number,  gas  flow 
angle,  static  pressure,  total  pressure  and  thickness  of  the  flow 
layer  of  the  gas  flow  behind  the  shock  wave  in  the  duct  respectively 
The  subscripts  1.  and  2.  respectively  represent  the  direction  of 
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gas  flow  in  the  front  and  back  of  the  cascade.  Using  the  flow 
parameters  in  the  front. and  behind  the  shock  wave  in  the  duct,  we 
can  write  the  momentum  conservation  equation  in  the  8  direction 
(through  continuous  equation  transformation) 

p  +  pwl  —  +  pWH/tQQtiji  —  ^  ^ 


Using  the  state  equation  and  the  sound  velocity  expression. 


we  get 


M  i  cos/?. 


—  —  O  +  —  rWMico»(0,  —  0) 


Equation  (2)  can  be  used  to  calculate  the  operating  attack 
angle  of  the  unstarted  choked  flow  pattern.  To  carry  out  the 
calculation,  it  is  also  necessary  to  be  supplemented  by  additional 
relations.  Let  us  assume  that  iii  can  be  used.  In  addition,  because 
the  region  behind  the  wave  is  the  upper  subsonic  region,  the  throat 
of  the  cascade  St  (or  the  effective  throat  created  due  to  viscosity) 
is  usually  in  a  critical  state.  Therefore,  the  following  equation 
can  be  used: 

—  u(  r  +  1’ 

M  \2  +(r  -  l)MV  (3) 

q(M)  is  defined  as  the  critical  contraction  ratio.  Under  the  con¬ 
dition  of  downstream  throttle,  the  throat  is  located  at  the  throttle. 

* 

In  calculation  the  extent  of  opening  of  the  throttle  A2  should  re¬ 
place  Sfe. 

Then  from  assumption  ii,  the  flow  direction  S  behind  the 
shock  wave  in  the  duct  is  the  average  geometric  direction  of  the 
cascade.  It  is  obtained  under  the  equi-fl  condition  that: 

Using  Equations  (2)-(4),  we  can  obtain  the  four  unknowns  B1(M1), 

IVKM^),  B(MX)  and  $2(M^)  with  given  Q  and  St/t  Cor  A^Ct)  where 
B^(M^)  is  the  operating  attack  angle  which  we  want  to  determine. 

During  the  process  p^/p2  is  set  ahead  of  time  and  it  is  checked 
against  other  parameters  of  the  cascade  when  the  solution  is  obtained. 
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Therefore,  it  is  an  iterative  process.  After  the  flow  parameters 
behind  and  in  front  of  the  shock  wave  in  the  duct,  the  static 
pressure  use  and  the  total  pressure  recovery  coefficient  of  the 
shock  wave  in  the  duct  can  be  directly  obtained  from  the  continuous 
equation 

P_  _  pT  _  b,w ,  cos;?,  —  f I  +  (r  ~~  QAf  ?/2l  t 

Pt  bwcotft  b  M  co$(l  L|  +  Ct  OMV2-*  \31 


Let  us  take  b  =  (b^+b2)/2  in  approximation  and  then  we  get 

p>  2 Q  fl  +  Cr  -  l)M72li£=n- 

+  \  Mcoafi  [1  +•  Cr  -  l)Mi/2]  ^-0  1 


For  the  gas  flow  behind  the  shock  wave  in  the  duct  in  the 
subsonic  cascade,  its  optimal  compression  efficiency  can  be  anal¬ 
yzed  using  the  results  of  dimensional  analysis  given  in  [10]. 

/  1  ~  n.y  _  L*-n 

\  0.14ij,  /  yi-Bjjt.2*  (y) 

nQ  is  the  varying  efficiency.  The  three  parameters  on  the  right 
hand  side  of  Equation  (7)  are 


L 


cos  /?,COSJ<£ 


Re 

oa-,*co*<^ 


]• 


U 


C8) 


where  a  =  c/t  is  the  density,  B  is  the  installation  angle,  <p  is 

s 

the  middle  arc  angle  and  Re  is  the  Reynolds  number  of  the  arc  length. 


In  order  to  consider  losses  due  to  shock  wave  in  the  duct  and 
boundary  layer  interference,  etc.,  a  correction  factor  f  is  intro¬ 
duced  to  the  relation  between  nQ  and  the  total  pressure  recovery 
coefficient : 


Reference  [7]  gathered  lots  of  data  and  obtained  the  follow¬ 
ing  using  a  linear  relationship 


/_l  11 

5  i  +  fjiM}  M,J 


(10) 


33 


x2  is  the  thickness  of  the  rear  fringe,  t/(t+i2)  can  be  used  to 
represent  the  expansion  area  ratio  whose  effect  on  the  loss  is 
reasonable.  Under  the  condition  that  a  profound  study  is  not 
carried  out,  we  are  going  to  temporarily  use  this  correction  fac¬ 
tor  expression  here. 


Using  Equations  (6)-(10),  we  can  write  the  total  pressure 
recovery  coefficient  of  the  supersonic  cascade 

pi  ri  +  Cr  -  l)/ij,WV2ir^Tf  2Q  M, co. A  /l  +  (r  - 

P?“  l~  +Cr-  DWV2  J  +  l  M  c«*  \l+  Cr  “  l)M{/2 )  \  C1 


Other  properties  of  the  cascade  are 


(12) 

Pi  pi  f  1  +  (r  —  i~)M!/2  ifSt 

P.  P?  L I  +<r-  OWj/2  J 

(13) 

£•  piM  -  0M?/2  irh 

*  P?  1  ‘  +  (r  ~  OMi/2  J 

(14) 

When  using  Equations  (13)  ,(1*0,  we  must  know  M2»  It  has  been 
known  that  [Ocotfijeotfi,]  —  I/>Wpiw<1  —  K  .  and 
Combining  the  two  equations,  we  can  get 

Cl  5) 

q(M)  is  called  the  flow  density  function  or  the  critical  contrac¬ 
tion  ratio,  which  is  defined  as 


*(M) 


A *  (  r  +  l  \£=fr 

A  "  M  V2  +  (r  -  OW1  J 


(16) 


Here  parameters  carrying  the  *  symbol  are  called  critical 
parameters . 

The  value  of  p^/Pi  obtained  from  Equation  (14)  should  be 
consistent  with  the  assumed  value*—  lotherwise ,  we' must  correct  the 
assumed  value  and  begin  the  iteration  process. 
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III.  SOLUTION  PROCEDURE  AND  RESULTS 
OP  CALCULATION 


It  is  not  complicated  to  solve  the  above  equations. 

Equations  (2)  -  (4)  may  be  used  individually  to  find  fa* 
and  correspondingly  to  obtain  the  gas  flow  parameters  behind 
the  channel  shock  wave.  In  obtaining  the  solution,  one  should 
first  assume  a  value  for  the  ratio  fiJpi  of  the  densities  in 
front  of  and  behind  the  blading.  It  is  best  in  actual  procedures 
of  solving  the  equation  to  assume  fa  first,  and  then  to  find 
the  root  of  from  Equation  (4).  There  are  two  roots,  cor¬ 
responding  to  the  supersonic  solution  and  the  subsonic  solution. 
Thus  there  is  one  fa  corresponding  to  two  M.  .  Equations  (5)- 
(16)  are  used  to  compute  the  other  properties  of  the  blading  and 
to  check  the  assumed  value  of  *>i/p»  .  This  is  an  iteration 

procedure.  When  the  necessary  accuracy  is  reached,  the  angle  of 
attack  and  other  properties  of  the  corresponding  supersonic  and 
trans-sonic  blading  may  then  be  obtained. 

The  results  obtained  by  using  this  method  on  the 
triangular  blading  in  Figure  3  are  shown  in  Figures  4-7.  In 
Figure  3,  the  shape  of  the  blading  and  necessary  data  are  given. 
This  is  also  the  blading  calculated  by  Fabri.  Since  our  model 
is  basically  one-dimensional,  its  results  may  be  borrowed  for 
blading  with  similar  geometric  and  dynamic  parameters. 

Since  experimental  data  are  lacking,  we  take  K  =  0.8 
from  Figure  2.  The  results  of  this  paper  are  computed  with  K=0.8. 
If  we  take  K  *  0.75,  the  curve  for  the  angle  of  attack  will  lie 
above  the  curve  for  K  =  0.8  (not  shown). 

The  calculated  results  for  3  types  of  one-dimensional  models 
of  Fabri  are  plotted  in  Figure  4.  Compared  to  results  of  our 
paper,  the  trends  of  the  variation  agree  and  the 


35 


Figure  3.  Parameters  ofl  the  triangular 
cascade  and  the  velocity  triangle. 

>— 32.4mm  r  -  60mm  S,  —  14.1mm 
«, -4*  0,-60* 


1 — average  geometric  direction 


Figure  4.  Calculation  of  the  operating  attack  angle  of  unstarted 
choked  regime  of  a  triangular  cascade  (K  =  0.8) 

1 — shock  wave  formed  in  the  average  direction  of  the  end  of  the 
duct  (.recommended  by  Chauvin);  2 — shock  wave  deflection  (recommend¬ 
ed  by  Fabri);  3 — formation  of  normal  shock  wave  in  the  direction  of 
the  incoming  flow;  4 — experimental  data  (ONERA  special  permission 
[4];  5 — calculated  results  in  this  work  (Q  ■  1.0) 
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difference  in  the  value  of  the  operating  attack  angle  is  not 
large.  Because  we  did  not  find  any  experimental  data  on  the  cal¬ 
culation  of  triangular  cascade,  an  experimental  curve  was  also 
included  as  a  reference  in  the  figure  which  was  allowed  to  be 
published  with  the  special  permission  of  ONERA,  Figure  5  re¬ 
flected  the  effect  of  the  axial  flow  density  ratio  Q'  on  the  oper¬ 
ating  attack  angle  of  the  unstarted  choked  cascade.  The  effect 
of  Q  has  not  been  considered  by  other  methods.  Q  >  1  corresponds 
to  the  increasing  contraction  ratio  at  the  throat  of  the  cascade. 
At  this  time,  it  requires  a  larger  variation  for  (3^  to  make  the 
starting  progress  along  the  choked  curve.  Figure  6  gave  the 
effect  of  throttle  behind  the  cascade  and  made  a  comparison  with 
the  result  obtained  from  a  two  dimensional  model  in  E5]-  The 
trend  of  variation  appears  to  be  consistent.  The  attack  angles 
only  differ  by  2-3°.  Part  of  the  reason  may  be  attributed  to  the 
fact  that  the  two  calculations  were  made  based  on  different  cas¬ 
cades.  The  two-dimensional  model  in  [5]  cannot  calculate  the  tran 
sonic  region  enar  =  1.  Therefore,  it  is  connected  by  a  dotted 
line  in  the  figure.  Figure  7  plotted  other  properties  of  the 
unstarted  choked  cascade,  including  the  total  pressure  loss  coeffi 
cient  U,  static  pressure  rise  P2/Pi  anc*  deflection  angle  A(3. 

Their  variation  trend  with  and  effects  on  Q  more  or  less  re¬ 
flected  the  actual  situation  [11,12]. 

The  calculation  of  loss  due  to  shock  wave  in  the  duct  indi¬ 
cated  that  the  total  pressure  loss  is  larger  than  the  total 
pressure  loss  caused  by  the  incoming  flow  M-^  and  the  normal  shock 
wave.  This  is  due  to  that  the  actual  wave  series  at  the  inlet  of 
the  cascade  should  contain  an  infinite  number  of  wavefronts  ahead 
of  the  quasi-normal  shock  wave  in  the  duct  produced  by  the  "pre¬ 
ceding  blades"* 
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FIGURE  5.  The  calculation  of  the  operating 
attack  angle  of  unstarted  choked  regime  of  super-  °Jf- 
sonic  and  transonic  compressor  cascades  (varying 
axial  velocity  density  ratio)  (K&0-.-8}  6>&)  • 
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LO  U  1.4  U  IS  2.0 

'FIGURE  7.  Calculation  of  properties 
of  the  unstarted  choked  regime  of 
supersonic  and  transonic  cascades 
(varying  axial  velocity  dens,  ratio) 


(FIGURE  6.  Calculation  of  unstarted  choked  regime 
jof  supersonic  and  transonic  planar  cascades  (var 
inq  throttle  ratio).  Calculated  results  in  [5]. 
Calculated  results  in  this  work. 


IV.  CONCLUSIONS 

This  paper  provided  a  relatively  simple  calculation  method 
which  is  capable  of  calculating  the  operating  attack  angle  and  the 
corresponding  characteristics  of  the  cascade  of  unstarted  choked 
regime  of  the  supersonic  and  transonic  cascades.  The  results  were 
compared  to  those  reported  in  C 33— C 5 ]  which  indicated  that  the 
method  introduced  in  this  paper  is  valuable.  It  can  also  consider 
the  effect  of  downstream  throttle  and  axial  velocity  density  ratio 
variation.  Q  has  a  very  important  effect  on  the  characteristics 
of  the  cascade  and  this  is  well  known.  More  calculation  should  be 
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carried  out  on  these  supersonic  and  transonic  cascades  with  some 
experimental  results  in  order  to  further  'improve  and  verify  this 
method.  It  is  more  important  to  investigate  the  selection  of 
the  value  of  K. 

During  the  derivation  and  calculation  stages  of  this  paper, 
instructions  have  been  given  by  Professor  Wu  Chung  Haw, 
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THE  NUMERICAL  COMPUTATION  FOR  DIFFRACTION  AND  REFLECTION  OF  A 
TRAVELING  SHOCK  WAVE  FROM  AN  OPEN  END  OF  A  SHOCK-TUBE 

Du  Xun,  Wei  Zbong-lei,  Li  Wen-xuan,  Wang  Jia-jun 

ABSTRACT 


When  a  shock  wave  travella  into  atmosphere  from  an  open  end  of  a  shock-tube, 
an  over- pressure  field  ia  induced  and  wavo  diagrama  of  shock  diffraction  and  reflec¬ 
tion  form  immediately.  In  general,  one  ean  always  compute  the  outflow  problem  by 
solving  the  unsteady,  two-dimensional  gatdynamic  equations.  Computing  this  problem 
numerically,  the  authors  take  advantage  of  Eulerian  coordinates,  fictitious  viscosity 
and  finite  difference  method.  In  order  to  decide  the  boundary  conditions  for  com¬ 
putation  at  the  end  of  tube,  first  we  find  out  analytical  solution  for  the  unsteady 
one-dimensional  flow  in  a  shock-tube,  sod  then  the  variations  of  the  velocity,  pres¬ 
sure  and  density  of  the  flqw  st  the  open-  end  of  the  tube  with  tima.  The  present 
numerical  results  are  in  agreement  with  the  experiments. 


The  wave  diagram  and  the  over-pressured  field  Induced  by 
traveling  shockwave  from  the  open  end  of  a  shock-tube  can  be  sim¬ 
plified  into  a  two-dimensional  diffraction  and  reflection  problem 
of  an  unsteady  shockwave.  This  set  of  gas  dynamic  equations  is 
non-linear.  The  solution  would  show  interruptions  with  time.  Even 
if  the  initial  value  is  a  continuous  function,  there  may  be 
interruptions  in  the  solution.  This  causes  some  difficulty.  In 
this  paper,  the  authors  adopted  Eulerian  coordinates,  artificial 
viscosity  and  a  finite  difference  method  to  obtain  the  solution. 

In  order  to  write  the  boundary  conditions  at  the  end  of  the  tube, 
an  analytical  solution  to  a  one-dimensional  gas  flow  in  a  tube  was 
first  found.  From  this  solution,  the  time  varying  curves  of  the 
velocity,  pressure  and  density  at  the  end  of  the  tube  can  be  com¬ 
puted.  The  calculated  results  were  compared  with  the  experimental 
ones,  and  the  two  were  found  to  be  consistent.  We  have  investigated 
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the  wave  forms  of  the  diffraction  of  the  shockwave  traveling  from 
the  end  of  the  tube,  the  slanted  reflection  of  the  shockwave  from 
the  ground,  the  normal  reflection  of  the  shockwave  from  a  perpen¬ 
dicular  barrier,  the  slope  climbing  process  of  the  shockwave  and 
the  formation  and  development  of  the  rotating  ring  at  the  end  of 
the  tube  using  a  shadow  tracing  method.  Simultaneously,  the  over¬ 
pressure  field  distribution  outside  the  end  of  the  tube  induced  by 
the  series  of  shockwaves  was  also  measured.  Numerical  computations 
were  also  carriedout  in  this  paper  for  the  diffraction  and  reflec¬ 
tion  flow  fields  after  the  shockwave  traveled  through  the  end  of 
the  tube. 


I.  Gas  dynamic  equations 


Neglecting  viscosity  and  thermal  conductivity  and  with  con¬ 
stant  specific  heat,  the  set  of  dimensionless  gas  dynamic  equations 
for  unsteady  gas  flow  is : 
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where  P^fi/p**  p-'p/p*,  *  —  5/««,  u  —  5/«„  * -—?/£),  y  —  y/D,  I— >i/(D/ J,), 

q*mm  specific  beat  of  air  *—1.4,  p,  and  pQ  are  the  den¬ 

sity  and  pressure  at  standard  atmospheric  conditions,  D  is  the 
characteristic  diameter  of  the  shock-tube.  The  symbols  with  a  bar 
on  top  represent  quantities  with  dimensions.  For  planar  problems 
v  ■  0;  for  symmetric  axial  problems,  v  ®  1. 
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As  shown  in  Figure  2,  analytical  solutions  can  all  be  found 
in  regions  (1)-C7)  for  the  gas  flow  in  a  one-dimensional  shock- 
tube.  At  the  end  of  the  tube,  the  flow  is  supersonic.  Therefore, 
the  velocity,  pressure,  density  at  the  end  of  the  tube  are  com¬ 
pletely  determined  by  the  one-dimensional  flow  inside  the  shock- 
tube.  Region  Cl)  is  an  undisturbed  region:  p1  ■  ■  1,  u^  ■  0. 

Region  (.2)  is  the  shockwave  compression  region: 


*-<70*  -l)/6 

(8> 

*,-5(0-O/0)l/« 

(9) 

p,  -  6070  +  00 

Oo) 

0  is  the  dimensionless  shockwave  moving  velocity.  It  is  defined 
as  0  =  0/aQ.  Region  (3)  is  the  uniform  region  behind  the  contact 
interruption  plane: 


ft  —  pi,  —  «i 

Pi  *■*  ?«/*( 

-  1  -  W5) 


(ID 

02) 

03) 


Here  «  —  «/«,  ;  therefore,  •  Region  (4)  is  the  high 

pressure  region  of  the  shock-tube.  In  the  computations,  let  p^  = 
Pjj  =  10,  and  Ujj  ■  0.  In  addition,  there  is  a  relation  between  p^ 
and  p2 : 


P*  ™  Pi  (l 


04) 


If  the  value  of  p^  is  given  for  a  flow  in  the  shock-tube,  the 
shockwave  intensity  p2  and  02  produced  can  be  calculated. 

In  regions  (5),  (6)  and  (7)  p/p1'2*  *  Ph-2^  *  constant. 

p  H 

Using  p  »  p  to  substitute  into  it,  we  get 

Cl 

p  —  p^t  P  —  05) 


43 


Therefore,  for  the  solutions  for  the  non-uniform  regions  C5), 
(6)  and  (7),  we  only  have  to- write  the  expressions and  u. 


Region  (5)  is  the  central  dilute  wave  region.  Its  solution 


is 


06) 


where  *“*//»  f/(//«,),  /  is  the  length  of  the  high  pressure 
section;  x  is  the  time  period  from  the  membrane  breaking  moment  of 
the  shock-tube. 


Region  (6)  is  the  non-uniform  dilute  wave  interference  region. 
Its  solution  is  (define  £  *  p/5) 


3  P  —  6(V  +  Of1  +  3V  +  2V  +  3 


07) 

08) 


Region  C7 )  is  the  simple  wave  region.  Its  solution  is 
«  -  (»/5)  -  «,  -  W5)  -  t  09) 

i.e.,  £  =  a-c  and  c  is  a  constant.  Substituting  Into  Equations 
Cl7)  and  (18),  we  obtain  the  equation  of  the  boundary  life  CF  of 
region  (.6)  and  region  (7). 


r  _  f(<)  _  4QV  -  IV  -  12r(V  -  1).  +  3(c*  -  1 )» 
,  -  JCi)  -  U  -  o  (5f<o  -  l£±±f±±=^] 


(20) 

(21) 


Because  the  first  family  characteristic  curve  is  a  straight 
line,  therefore 

«  —  *(«)•■(»  + «Xr  —  *(«))  (22) 


The  solution  for  region  (.7)  can  be  obtained  from  the  simul¬ 
taneous  Equations  (19)-(22). 
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Prom  Equation  (16)  we  get  £  *  u/5  =  1-a.  Substituting  it 

into  Equations  (17)  and  (18),  we  obtain  the  equation  x  =  1  +  5t  - 
2/? 

6x  J  for  the  boundary  line  BC  between  region  (6)  and  region  (5). 
The  coordinates  of  point  C  can  also  be  obtained 

r.--,  x.-1  +  i-A  (23) 

«l  «!  . 

The  boundary  line  CD  between  region  (3)  and  region  C7)  is  a 
characteristic  line.  Its  equation  is 

*  —  x,  “  («, +  «»)(*  —  O  “■  («» +  ~  *«)  •  (24) 

The  equation  of  the  contact  interruption  plane  AD  is 

*  -  1  +  «, r  (25) 

Solving  Equations  (24)  and  (25)  simultaneously,  we  get  the 
coordinates  for  point  D: 

vo  “  "  l  +  ^  (26) 

*\  a. 


In  the  case  that  it  is  supersonic  at  the  end  of  the  tube  and 
the  length  of  the  shock-tube  x,  —  L/K*0  Cl,  L  are  the  length  of 
the  high  pressure  section  and  the  shock-tube  respectively),  for 
given  Initial  values  of  regions  (.1)  and  (JO,  the  entire  flow  field 
in  the  shock-tube  can  be  obtained  by  using  the  equations  for 
regions  (2),  (3),  (5)  and  (.7)  as  described  above.  Our  experimental 
condition  satisfied  the  conditions  that  the  flow  was  supersonic  at 
the  end  of  the  tube  and  */<r0*  Therefore,  the  above  equations 
can  be  used  to  calculate  the  distribution  of  each  parameter  when 
the  shockwave  reaches  the  end  of  the  tube  (.Figure  3(a)).  Then 
the  variations  of  the  flow  parameters  at  the  end  of  the  tube  were 
converted  from  the  time  beginning  with  the  breaking  of  the  membrane 
r—f/  O/a0)  i  to  the  time  /  —  i/CD/ 'aQ ) .  Prom  the  figure  it  can  be 

seen  that,  if  we  let  then  the  calculated  time  of  the 

arrival  of  the^shgckwave  at  the  end  of  the  tube  *  4.065,  * 

5.511,  x 2  *  0o5~,  and  *  7.5.  Based  on  this,  we  can  also  deduce 


■<T 
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the  relation  between  the  dimensionless  time  t  and  x  to  be 


/  -  (r  -  rt)<//D)  -  (r  -  4.065)  •  2 


The  practical  procedure  used  in  this  paper  to  determine  the 
boundary  conditions  at  the  end  of  the  tube  for  numerical  computa 
tions  is:  Starting  at  t  =  0 ,  the  variation  of  each  parameter  at 


the  end 

of  the 

tube  with  time  is: 

<* 

1) 

when 

0  <  t  <  4, 

,  it  belongs  to  region  (2): 

ft-  10,  -2.1775,  p, -3.81 

2) 

when 

/,</</, 

,  it  belongs  to  region  (3): 

ft— 10,  ft—  2.1775,  ft— 31.38,  —  0.5646 

> 

i 

3) 

when 

h  <  t  <  |, 

u  = 

,  it  belongs  to  region  (7): 

5(a-0,129) 

The  variation  of  a(t)  is  obtained  from  Equation  (22)  and  by 

noticing  x  *  (t+8.13)/2: 

I  -  2  [f(«)  +  -  8.13  (27) 

l  6  •  —  0.6454 

where  the  expressions  for  ~(a)  and  x(.a)  are  obtained  from  the 
equation  of  CF  after  substituting  the  value  <r  — W5) —  0.129 
into  Equations  (.20)  and  (.21),  i.e., 
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(20 

(22') 


r 


r  -  *(„)  -  (0.3624  +  0.1903«  -  0.475«*)/«* 

N  ,  «  ,,a*  f «,/  v  «*  +  0.129«  +  0.4917  1 

x  ■“  *(j)  “(i  -  0.129)  ^5f(«)  —  j;  J 

For  convenience  in  using  the  computer,  after  calculating 
the  numerical  relationship  between  a  and  t  for  Equation  (.271, 
it  is  not  difficult  to  obtain  an  approximate  equation  as  follows: 

*  —  0.6S63  -  0.017661,  m  -  2.6363  -  0.08831  (28) 

The  density  p  and  pressure  p  are  calculated  using  Equations 
(15)  and  (16). 

4 )  When  t  >  t3,  it  belongs  to  region  (6).  Using  x  =  13 
and  t  =  (t+8.13)/2  to  substitute  into  Equations  (20)  and  (21), 

we  find  that  the  variations  of  u  and  a  with  time  are  as  follows: 

/  -  [35«  -  6(«*  +  l)*1  +  U*  +  2«*  +  3l/(4«»)  -  8.13  (29) 

13  -  f[15r  -  (26«»  +  30)5*  +  3«4  +  6 j*  +  15]/(8«0  (30) 

From  Equation  (.30)  £  and  a  can  be  solved.  By  substituing 
them  into  Equation  (29),  the  corresponding  t  value  can  be  computed. 
After  performing  the  numerical  calculations,  the  approximate 
equation  can  then  be  found  for  convenience  in  using  the  computer. 

In  order  to  further  improve  the  approximation,  the  approximate 


equation  is  chosen  by  sections  as 

follows : 

«  —  0.5341  +  0.0015881  -  0.0002123/*  1 
m  —  3.285  —  0.2248/  +  0.00611/*  J 

a  6.87  </<  12.7  It 

(31) 

4  —  0.5709  —  0.004425/  4*  0.00003313/* 
m  —  2.781  —  0.1423/  +  0.002739/* 

}  a  12.7  <  /  <  18.61  at 

(32) 

•  —  0.5726  -  0.004573/  +  0.00003604/* 
u  —  2.1165  -  0.07136/  +  0.0008482/* 

J  a  18-61  </<  33.43  It 

(33) 

•  —  0.5528  —  0.003405/  4-  0.00001884/* 
it  -  1.433  -  0.02933/  +  0.000203/* 

|  a  33.43  </<  57.31  It 

(34) 
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After  a  and  u  are  solved,  then  the  corresponding  p  and  p  are 
derived  from  Equation  (15).  The  variations  of  various  parameters 
at  the  end  of  the  tube  with  time  are  shown  in  Figure  3(b). 

It  should  be  noticed  that  when  t  =  57.31,^  =  5*  S  *  0.4195 
=  a.  Therefore,  when  t  >  57.31*  a  subsonic  flow  begins  to  appear 
at  the  end  of  the  tube.  The  outside  of  the  tube  would  influence 
the  inside,  and  the  solution  for  region  (6)  becomes  ^inapplicable . 
Henceforth,  the  computation  carried  out  in  this  paper  will  not 
exceed  a  maximum  time  period  t  =  57-31. 

3.  Initial  conditions 


When  t  =  0,  with  the  exceptions  that  —  10,  v,  —  o,  *,  —  2.1775, 
2-;  and  p2  *  3.81  at  the  end  of  the  tube,  in  the  remaining 
regions  we  will  assume  p  *  p  =  1,  and  *«*»»  —  0, 


III.  Difference  equations 

1.  The  format  of  difference  for  the  internal  points  of  a 
planar  problem  in  gas  dynamic  equations  [ 2 J 

A  rectangular  lattice  is  formed  in  a  space  coordinate  system 
(.x,y,t)  with  unit  length  Ax  =  h1.  Ay  *  b2  and  At  *  x.  Let  the 
parameter  A  at  a  lattice  point  with  coordinates  Cmh1,  lh2>  nx) 
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be  expressed  as  A^  Furthermore,  let  us  use  the  abbreviation 
A  for  it.  In  addition,  let  us  decide  to  ignore  all  the  other 
parameters  with  the  same  indices  m,n,l.  The  set  of  gas  dynamic 
equations  after  talcing  the  fictitions  viscosity  into  considera¬ 
tion  can  be  transformed  into  the  following  set  of  difference 
equations: 

r**  -  r  -  4*  -  f— >  -  41  \  i -  /.) 

2  Z  2 

-  o.  -  /—,)]  +  j  -  /,)  -  -  /,-,)]  05) 

where  the  definitions  of  f,F  and  G  are  shown  in  Equations  (2)- 
(4)  and 

*•*1  ™  ~  (°»*i  +  ®»)»  At)  "  ~  (Ati  +  ft)  (36) 

Furthermore,  Ai  —  r/A,,  ^  —  r/A„  Aj/A,  —  tg<£  a  and  8  are  defined  as: 

a  —  (offsin1^,  0  —  uffcosty  (37) 

where  »-*(«  +  -),  ^-V^T+Ti,  *-(«*+ •-VpTp  a called 

the  Courant  number  at  the  nodal  point  of  the  lattice,  ;!r  is  the 
stabilization  constant  and  the  stabilization  condition  of  the 
difference  equation  requires  [2]: 

«r l/<r  (38) 

Since  the  stabilization  condition  (38)  is  valid  at  any  lat¬ 
tice  point  under  any  condition,  we  can  let  crj ■“  max  <*!.<.  .  if  <r;«£i 
then  the  following  equation  must  also  be  valid: 

o;<u<i/9i  (39) 

In  the  numerical  computation  of  the  computer,  it  can  be  deter¬ 
mined  that  each  level  has  the  same  a  0  value.  Based  on  the  equation 

k  —  »i/(f  +  •>«« 

the  corresponding  value  of  k  at  each  level  can  be  obtained.  Since 
the  values  of  h^  and  h2  have  been  selected,  the  corresponding  x 
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value  at  this  level  can  then  be  determined.  Thus  x  is  selected. 

6) 

Any  ^  which  satisfied  *;<»<  l/*f  must  also  satisfy  the  condi¬ 
tion  l/«f  at  any  nodal  point  in  the  lattice  at  the  same 

level.  In  the  actual  numerical  computation,  the  selected  value 

of  at  every  level  is  the  same  constant  value  and  the  selected 
° 

value  of  j#' is  also  the  same  for  all  levels.  But  at  every  nodal 
point  of  the  lattice  at  every  level  the  value  of  a  is  different 
which  is  determined  by  the  equation  »  —  .  The  choice  of  t 

in  the  machine  is  automatic. 

2.  The  format  of  difference  at  the  boundary  point  in  a 
planar  problem 

From  the  difference  format  of  internal  point  discussed  in  the 
previous  section,  it  is  understood  that  the  calculation  procedure 
uses  the  known  parameters  at  five  points  (M,l),  (M,i  ±  1)  and 
Cm  ±  1,  l)  at  nth  level  to  compute  the  unknown  parameters  (u,v, 
p,p)  at  point  (m,l)  of  level  n  +  1.  Equation  (35)  does  not 
apply  to  boundary  points.  Some  modification  is  necessary  with 
regard  to  the  difference  format  of  the  boundary  points. 

For  horizontal  solid  wall  (including  horizontal  symmetric 
axis)  and  vertical  solid  wall,  we  adopted  the  method  which  trans¬ 
fers  the  boundary  difference  format  into  the  difference  format  for 
internal  points.  The  actual  procedure  is:  For  horizontal  solid 
wall,  assuming  that  the  index  number  is  1  in  the  y-direction  of 
the  solid  wall  and  then  sending  to  the  unit  of  and  then 

changing  sign.  For  the  other  parameters  »/><*•»  •»«*•  ,  their 

original  values  are  transferred  to  p»*i»  *»»*.  *t*i  respectively. 

This  method  in  essence  is  to  reflect  the  value  of  the  /  ±  1  row 
into  mirror  image  to  the  / T 1  row.  After  a  row  of  points  is 
generated,  the  original  horizontal  solid  wall  boundary  can  be  con¬ 
sidered  as  internal  points.  Computation  can  be  carried  out  using 
the  internal  point  difference  format  Equation  (35).  Similarly, 
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mirror  reflection  is  made  for  the  vertical  solid  wall  which  is 
to  transfer  to  with  a  sign  change  -and  to  transfer 
the  original  values  of  the  other  parameters. 


3.  Difference  format  of  a  slant  wall 
Generally 

9e»«aliy,  the  difference  format  of  a  slant  wall  will  attempt 
to  make  the  slant  wall  to  coincide  with  the  diagonal  line  of  the 
difference  lattice.  Thus,  we  can  rotate  the  slant  wall  by  an 
angle  to  transform  the  slant  wall  into  a  "horizontal"  solid 
wall.  After  the  difference  format  of  this  "horizontal"  solid  wall 
is  written,  we  can  then  transfer  the  coordinate  system  to  obtain 
the  difference  format  of  the  slant  wall: 


1'*'  “■  ?" - ^  CO«  (+1  *”  Pm- 1.1-*)  (eottydmu+i  +  datyGm-u) 

2  co*  <p 

J+  (40 


l\  \ 

(1  ] 

/  %  1 

tm  \ 

1 

+ 

w 

H 

it’**, 

\p2*  +  5pm  / 

j.  3  —  j 

| 

W  +  V>2/ 

Symbols  with  the  sign  indicate  that  they  are  the  para¬ 
meters  in  the  (x,y)  coordinate  system  which  is  formed  by  rotating 
the  original  coordinate  system  by  an  angle  ip  with  x  axis  parallel 
to  the  slant  wall  and  y  axis  perpendicular  to  the  slant  wall. 
Their  relations  with  those  in  the  original  coordinate  system 
(x,y)  are: 


2  mco*</»  4- vsin^,  5  ™  »co*<> 


(43) 


On  the  slant  wall  ?  —  0.  Substituting  this  into  the  above 
equation,  we  get 

u  *  u cotip,  »  —  3»in  <p  (44) 
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Based  on  the  above  equations ,  the  values  of  parameters 
(>, «»  and  •»  on  the  slant  wall  can- finally  be  obtained. 

M .  The  difference  format  in  axial  symmetry  problems 

The  coordinate  x  represents  the  direction  of  the  symmetric 
axis  and  the  coordinate  y  axis  represents  the  direction  perpendi¬ 
cular  to  the  symmetric  axis.  After  adding  another  term  to 

the  planar  difference  format,  its  internal  point  difference  format 

iD  P*1  -  T  -  (F.*  -  F_)  -  (Cm  -  G,_)  -  rV  . 

2  2 

“t  y  /«-f)  *“  “*  fm — ill 

+  J  -  //)  -  -  /,-,)]  (45) 

where  the  definitions  of  f,  P,  G  and  V  are  shown  in  Equations  (2)- 
(5).  However,  on  the  symmetric  axis,  because  y-*-0  and  v-K),  the 
difference  format  is 

r*  p—)  -  2wln  +  i.  [a.+i<j.+1  -  /.) 

2  2 

/•— <)  1  +  pl+iQl-H  ~~  it)  (46) 


IV.  Analysis  of  the  computed  results 

We  have  computed  four  different  conditions:  Axial  symmetric 
flow  with  vertical  barrier  wall  without  ground;  two-dimensional 
planar  flow  with  slant  wall  of  actual  dimensions  with  ground;  and 
two-dimensional  planar  flow  within  two  small  areas  with  ground  and 
slant  wall.  The  objectives  of  the  first  three  computations  are  to 
find  the  correction  coefficients  of  the  three-dimensional  effect 
and  the  ground  reflection  effect.  The  last  computation  was  done  to 
examine  the  influence  of  the  choice  of  the  lattice  boundary  on  the 
calculated  result. 
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In  computation,  the  values  of  the  top  and  bottom  boundary 
points  at  the  end  of  the  tube  should  be  The  varia¬ 

tions  of  p  and  p  values  with  time  should  be  the  same  as  the  values 
of  the  points  inside  the  tube.  In  the  meantime,  please  notice 
that  inside  the  computed  area  it  is  a  free  boundary  which  is  that 
before  the  shockwave  propagates  through  p  and  p  and  1  and  u  and 
•  are  0  at  the  boundary.  After  the  shockwave  arrived,  the  follow¬ 
ing  simple  extrapolation  method  is  adopted:  Transferring  all  the 
computed  parameter  values  of  row  1  at  level  n  to  the  positions  of 
the  corresponding  points  in  row  l+l .  The  value  of  free  boundary 
in  row  l  at  level  n+1  can  thus  be  computed.  Here  l  is  the  maximum 
number  of  lattice  point  in  the  computation  in  the  y  direction. 

Based  on  the  above  computation  method,  some  of  the  results 
are  shown  in  Figures  4-6.  The  pressure  distribution  of  the  axial 
symmetric  flow  with  vertical  barrier  wall  is  shown  in  Figure  4. 
After  the  shockwave  has  traveled  through  the  end  of  the  tube,  the 
diffracted  shockwave  propagating  forward  will  create  a  normal 
reflection  when  it  collides  with  the  vertical  barrier  wall.  This 
normally  reflected  shockwave  will  also  be  compressed  by  the  high 
pressure  gas  flow  from  the  end  of  the  tube.  These  phenomena  are 
consistent  with  the  experimental  results  [1]. 


Figure  4.  Axial  syrrmetric 
.  condition 

(f„  m  lukg/cm',  t  ■  4.J2,  a  —  150) 


Figure  5.  Plane  condition 

(f„  -  lOkg/cm',  t  -  4.43,  a  -  150) 


The  pressure  distribution  of  the  planar  flow  with  a  vertical 
barrier  wall  is  shown  in  Figure  5.  Comparing  the  shockwave  diff¬ 
raction  pattern  of  the  axial  symmetry  flow  (Figure  4)  and  planar 
flow  (Figure  5)  at  the  same  time,  it  can  be  known  that  the  shock- 
wave  velocity  of  the  planar  flow  is  faster  than  that  of  the  axial 
symmetric  flow.  This  explains  that  the  induced  pressure  field  by 
the  axial  symmetric  flow  is  less  than  that  induced  by  the  planar 
flow. 


The  computed  result  of  the  shockwave  climbing  a  45°  tilted 
wall  (Figure  6)  and  the  observed  wave  series  diagram  [1]  are  also 
consistent , 


Based  on  the  calculated  data  using  the  axial  symmetry  and 
vertical  barrier  wall  method,  we  are  drawing  the  variations  of 
the  over-pressure  with  the  distance  from  the  end  of  the  tube 
radially  at  90°  and  150°  angles  (Figure  7).  They  are  also  more 
or  less  consistent  with  the  results  of  the  measurement  of  the 
over-pressure . 

<u 

£  02 
01 


I 

Figure  7 
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The  pressure  variations  with  time  are  drawn  (Figure  8)  for 
three  fixed  points  which  are  6m  high  from  the  ground  at  45°,  60° 
and  90°  angles  in  the  radial  direction  for  the  planar  computation 
condition.  The  effect  of  the  slant  wall  on  the  reflected  wave 
almost  cannot  be  noticed.  We  have  made  measurements  in  the  model 
experiment  at  90°  in  the  radial  direction  without  discovering  any 
effect  of  the  slant  wall  on  the  reflected  wave.  This  is  consistent 
with  the  computed  result. 

We  are  most  interested  in  the  diffraction  pressure  on  the 
rear  wave  at  150°C  in  the  radial  direction.  The  actual  computed 
result  after  taking  the  correction  coefficient  of  the  three-dimen¬ 
sional  effect  into  consideration  indicated  that  the  over-pressure 

p 

of  air  at  the  end  of  the  tube  is  10  kg/cm  and  the  free  diffrac¬ 
tion  over-pressure  at  the  corresponding  point  is  0.06.  The  model 
experiment  and  actual  on  the  spot  experiment  values  are  approxi¬ 
mately  0.0645.  This  again  explains  that  the  computed  value  is 
consistent  with  the  experimental  value. 

From  the  pressure  distribution  diagrams  4-6,  we  noted  the 
existence  of  a  low  pressure  nucleus  region  and  the  distribution 
of  equi-pressure  lines  surrounding  this  nucleus  at  the  upper  side 
out  of  the  end  of  the  tube.  It  is  equivalent  to  a  whirlpool.  This 
is  also  similar  to  the  rotating  ring  observed  experimentally  [1] 
although  the  position  of  the  two  is  slightly  different. 
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THE  NUMERICAL  COMPUTATION  FOR  DIFFRACTION 
AND  REFLECTION  OF  A  TRAVELLING  SHOCK  WAVE 

FROM  AN  OPEN  END  OF  A  SHOCK-TUBE 

* 

Da  Xun  Wei  Zhong-lei  Li  Wen-xuan  Wang  Jia-jun 

(Btijing  UiuvtTnti) 

Abstract 

When  a  shock  wave  travella  into  atmosphere  from  an  open  end  of  a  shock-tube, 
an  over-pressure  field  is  induced  and  wave  diagrams  of  shock  diffraction  and  reflec¬ 
tion  fora  immediately.  In  general,  one  can  always  compute  the  outflow  problem  by 
solving  the  unsteady  two-dimensional  gaadynamic  equations.  Computing  this  problem 
numerically,  the  authors  take  advantage  of  Eulerian  coordinates,  fictitious  viscosity 
and  finite  difference  method.  In  order  to  decide  the  boundary  conditions  for  com¬ 
putation  at  the  end  of  tube,  first  we  find  out  analytical  solution  for  the  unsteady 
one-dimensional  flow  in  a  shock-tube,  and  then  the  variations  of  the  velocity,  pres¬ 
sure  and  density  of  the  flow  at  the  open  end  of  the  tube  with  time.  The  present 
numerical  results  are  in  agreement  with  the  experiments. 


